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Abstract: 
O ' 

{Sj ■ We discuss in this Chapter a series of theoretical developments which motivate the 

introduction of a quantum evolution equation for which the eikonal approximation results 
in the geodesies of a four dimensional manifold. This geodesic motion can be put into 
correpondence with general relativity. The well-known problem of the self-interaction of a 
relativistic charged particle is studied from the point of view of a manifestly Lorentz co- 
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t — • ■ variant classical theory admitting the particle mass as a dynamical variable, i.e., the theory 
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is intrinsically off-shell. The evolution parameter, r, is an invariant parameter that can be 
identified with Newton's time (sometimes called, as in Schwinger's formulation, "proper 
q ' time" ) . Gauge invariance requires the definition of five gauge fields; the fifth field has for its 
source the matter density. Together with the results of Gupta and Padmanabhan, showing 
a connection between the radiation reaction force and geometry, this structure motivates 
an investigation into the connection between these dynamical equations and gravity. We 
show that the fifth gauge field can indeed be absorbed into a conformal metric in the kine- 
matic terms, which then results in a geodesic equation generated by the conformal metric 
and the standard Lorentz force. We then go on to show that the generalized radiation field 
passing through an optical medium with non-trivial dielectric tensor results in an analog 
gravity for the eikonal approximation for an arbitrary metric. A mathematically simpler 
system for which the eikonal approximation provides the geodesic motion on a four dimen- 
sional pseudo-Riemannian manifold is that of the Stueckelberg-Schrodinger equation with 
a spacetime dependent tensor (of the form of the Einstein metric tensor), somewhat 
analogously to a gauge field, coupling to the kinetic terms. This theory can be realized as 
a quantum theory in a flat spacetime, obeying the rules of the standard quantum theory 
in Lorentz covariant form. Since the geodesies predicted by the eikonal approximation, 
with appropriate choice of g^, can be those of general relativity, this theory provides a 
quantum theory which underlies classical gravitation, and coincides with it in this classical 
ray approximation. This result is the principal content of this work. 

In order to understand the possible origin of the structure of this form of the 
Stueckelberg-Schrodinger equation, we appeal to the approach of Nelson in constructing 
a Schrodinger equation from the properties of Brownian motion. Extending the notion of 
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Brownian motion to spacetime in a covariant way, we show that such an equation follows 
from correlations between spacetime dimensions in the stochastic process. 
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I. Introduction 

We shall deal with the problem of constructing a quantum evolution equation which 
describes, in the eikonal approximation, the geodesies of general relativity, on essentially 
three levels which include the basic motivations arising from the self-interaction problem of 
a relativistic charged particle in the framework of a generalized electromagetic theory, an 
analysis of the conformal structure that is associated with this theory, and the discussion of 
the eikonal approximation of wave equations that lead to our basic result. To understand 
the structure of the realtivistically covariant quantum equation, we study a relativistic 
generalization of Brownian motion, and follow Nelson's construction of the Schrodinger 
equation [1] for which local spacetime correlations result in a modification of the kinematic 
terms by introducing coupling to a Lorentz tensor g^ u somewhat analogous to a gauge field. 

It is quite remarkable that Gupta and Padmanabhan [2] , using essentially geometrical 
arguments (solving the static problem in the frame of the accelerating particle with a curved 
background metric), have shown that the description of the motion of an accelerating 
charged particle must include the radiation terms of the Abraham-Lorentz-Dirac equation 
[3], 

fj^ qnft A rp M 

where m is the electron mass, including electromagnetic correction, s is the proper time 
along the trajectory in spacetime, F£ is the covariant form of the electromagnetic 

force tensor, e is the electron charge, and 

2e 2 ,d 3 ^ d 2 x v d 2 x„dx» . 

3 c 3 [ ds 3 ds 2 ds 2 ds' [ ' ' 

Here, the indices /(/, z/, running over 0,1,2,3, label the spacetime variables that rep- 
resent the action of the Lorentz group; the index raising and lowering Lorentz invariant 
tensor n^v is of the form diag(—l, +1, +1, +1). The expression for T M was originally found 
by Abraham in 1905 [3], shortly after the discovery of special relativity, and is known as 
the Abraham four-vector of radiation reaction. Dirac's derivation [3] was based on a direct 
application of the Green's functions for the Maxwell fields, obtaining the form (1.1). 

Recognizing that the electron's acceleration precludes the use of a sequence of "in- 
stantaneous" inertial frames to describe the action of the forces on the electron [4], they 
carry out a Fermi- Walker transformation [5] , going to an accelerating frame (assuming con- 
stant acceleration) in which the electron is actually inertial, and there solve the Coulomb 
problem in the curved coordinates provided by the Fermi- Walker transformation. Trans- 
forming back to laboratory coordinates, they find the Abraham-Lorentz-Dirac equation 
without the direct use of the Maxwell Green's functions for the radiation field. This result, 
suggesting the relevance of curvature in the spacetime manifold, such as that generated 
by sources in general relativity, along with other, more elementary manifestations of mass 
renormalization (such as the contribution to the mass due to electromagnetic interactions 
and the identification of the Green's function singularity contribution with part of the elec- 
tron mass) , carries an implication that the electron mass may play an important dynamical 
role. 
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Stueckelberg, in 1941 [6], proposed a manifestly covariant form of classical and quan- 
tum mechanics in which space and time become dynamical observables. They are therefore 
represented in quantum theory by operators on a Hilbert space on square integrable func- 
tions in space and time. The dynamical development of the state is controlled by an 
invariant parameter r, which one might call the world time, coinciding with the time on 
the (on mass shell) freely falling clocks of general relativity. Stueckelberg [6] started his 
analysis by considering a classical world-line, and argued that under the action of forces, 
the world line would not be straight, and in fact could be curved back in time. He identi- 
fied the branch of the curve running backward in time with the antiparticle, a view taken 
also by Feynman in his perturbative formulation of quantum electrodynamics in 1948 [6]. 
Realizing that such a curve could not be parametrized by t (for some values of t there 
are two values of the space variables), Stueckelberg introduced the parameter r along the 
trajectory. 

This parameter is not necessarily identical to proper time, even for inertial motion for 
which proper time is a meaningful concept. Stueckelberg postulated the existence of an 
invariant "Hamiltonian" K, which would generate Hamilton equations for the canonical 
variables and p^ of the form 

*■ = (1.3) 



and 



where the dot indicates differentiation with respect to r. Taking, for example, the model 

*°-w (L5) 

we see that the Hamilton equations imply that 

M 

It then follows that 



^ = T7 ( L6 ) 



(1.7) 



dx. p 

~dt ~ E 1 

where p° = E, where we set the velocity of light c = 1; this is the correct definition for the 
velocity of a free relativistic particle. It follows, moreover, that 

V 1 = -*?r- = -m* 

With our choice of metric, dx^dx^ = —ds 2 , and p^p^ = — m 2 , where m is the classical 
experimentally measured mass of the particle (at a given instant of r). We see from this 
that 

^f! = at ( i9) 

dr 2 M 2 ' { } 
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and hence the proper time is not identical to the evolution parameter r. In the case that 
m 2 = M 2 , it follows that ds = dr, and we say that the particle is "on shell". 

For example, in the case of an external potential V(x), where we write x = x^, the 
Hamiltonian becomes 

K = ^M +V{X) (L10) 

so that, since K is a constant of the motion, m 2 varies from point to point with the 
variations of V(x). It is important to recognize from this discussion that the observable 
particle mass depends on the state of the system (in the quantum theory, the expectation 
value of the operator provides the expected value of the mass squared) . 

One may see, alternatively, that phenomenologically the mass of a nucleon, such as 
the neutron, clearly depends on the state of the system. The free neutron is not stable, 
but decays spontaneously into a proton, electron and antineutrino, since it is heavier than 
the proton. However, bound in a nucleus, it may be stable (in the nucleus, the proton may 
decay into neutron, positron and neutrino, since the proton may be sufficiently heavier 
than the neutron). The mass of the bound electron (in interaction with the electromag- 
netic field), as computed in quantum electrodynamics, is different from that of the free 
electron, and the difference contributes to the Lamb shift. This implies that, if one wishes 
to construct a covariant quantum theory, the variables E (energy) and p should be inde- 
pendent, and not constrained by the relation E 2 = p 2 + m 2 , where m is a fixed constant. 
This relation implies, moreover, that m 2 is a dynamical variable. It then follows, quantum 
mechanically, through the Fourier relation between the energy-momentum representation 
of a wave function and the spacetime representation, that the variable t, along with the 
variable x is a dynamical variable. Classically, t and x are recognized as variables of the 
phase space through the Hamilton equations. 

Since, in nature, particles appear with fairly sharp mass values (not necessarily with 
zero spread), we may assume the existence of some mechanism which will drive the parti- 
cle's mass back to its original mass-shell value (after the source responsible for the mass 
change ceases to act) so that the particle's mass shell is defined. We shall not take such 
a mechanism into account explicitly here in developing the dynamical equations. We shall 
assume that this mechanism is working, and it is a relatively smooth function (for example, 
a minimum in free energy which is broad enough for our off-shell driving force to work 
fairly freely) 1 . 

In an application of statistical mechanics to this theory [8], it has been found that a 
high temperature phase transition can be responsible for the restriction of the particle's 
mass (on the average, in equilibrium). In the classical theory, the non-linear equations 
induced by radiation reaction may have a similar effect [9]. 

1 A relativistic Lee model has been worked out which describes a physical mass shell 
as a resonance, and therefore a stability point on the spectrum [7], but at this point it is 
not clear to us how this mechanism works in general. It has been suggested by T. Jordan 
[personal communication] that the definition of the physical mass shell could follow from 
the interaction of the particle with fields (a type of "self-interaction"); this mechanism 
could provide for perhaps more than one mass state for a particle, such as the electron and 
muon and the various types of neutrinos, but no detailed model has been so far studied. 
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The Stueckelberg formulation implies the existence of a fifth "electromagnetic" poten- 
tial, through the requirement of gauge invariance, and there is a generalized Lorentz force 
which contains a term that drives the particle off-shell, whereas the terms corresponding to 
the electric and magnetic parts of the usual Maxwell fields do not (for the nonrelativistic 
case, the electric field may change the energy of a charged particle, but not the magnetic 
field; the electromagnetic field tensor in our case is analogous to the usual magnetic field, 
and the new field strengths, derived from the r dependence of the fields and the addi- 
tional gauge field, are analogous to the usual electric field, as we shall see). The second 
quantization of this gauge theory has been carried put as well [10]. 

In the following, we give the structure of the field equations, and show that the 
standard Maxwell theory is properly contained in this more general framework. Applying 
the Green's functions to the current source provided by the relativistic particle, and the 
generalized Lorentz force, we obtain equations of motion for the relativistic particle which 
is, in general off-shell. As in Dirac's result, these equations are of third order in the 
evolution parameter, and therefore are highly unstable. However, the equations are very 
nonlinear, and give rise to chaotic behavior [9]. 

Our results exhibit what appears to be a strange attractor in the phase space of 
the autonomous equation for the off-mass shell deviation. This attractor may stabilize 
the electron's mass in some neighborhood. We conjecture that it stabilizes the orbits 
macroscopically as well, but a detailed analysis awaits the application of more powerful 
computing facilities and procedures. 

We then show that the fifth (scalar) field can be eliminated through the introduction 
of a conformal metric on the spacetime manifold [11]. The geodesic equation associated 
with this metric coincides with the Lorentz force, and is therefore dynamically equiva- 
lent. Since the generalized Maxwell equations for the five dimensional fields provide an 
equation relating the fifth field with the spacetime density of events, one can derive the 
spacetime event density associated with the Friedmann-Robertson- Walker solution of the 
Einstein equations. The resulting density, in the conformal coordinate space, is isotropic 
and homogeneous, decreasing as the square of the Robertson- Walker scale factor. Using 
the Einstein equations, one see that both for the static and matter dominated models, 
the conformal time slice in which the events which generate the world lines are contained 
becomes progressively thinner as the inverse square of the scale factor, establishing a 
simple correspondence between the configurations predicted by the underlying Friedmann- 
Robertson- Walker dynamical model and the configurations in the conformal coordinates. 

The conformal metric is not, however, even locally equivalent to a Schwarzschild 
metric. To achieve a more general framework for achieving an underlying model for gravity, 
we study the eikonal approximation of the (generalized) electromagnetic equations in a 
medium with non-trivial dielectric tensor. 

It has been known for many years that the Hamilton-Jacobi equation of classical 
mechanics defines a function which appears to be the eikonal of a wave equation, and 
therefore that classical mechanics appears to be a ray approximation to some wave theory 
[12]. The propagation of rays of waves in inhomogeneous media appears, from this point 
of view (as a result of the application of Fermat's principle), to correspond to geodesic 
motion in a metric derived from the properties of the medium [13]. This geometrical 
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interpretation has been exploited recently by several authors to construct models which 
exhibit three dimensional analogs of general relativity by studying the wave equations of 
light in an inhomogeneous medium [14], and, to achieve four dimensional analogs, sound 
waves and electromagnetic propagation in inhomogeneously moving materials [15]. Visser 
et al [16] have pointed out that condensed matter systems such as acoustics in flowing 
fluids, light in moving dielectrics, and quasiparticles in a moving superfluid can be used to 
mimic kinematical aspects of general relativity. Leonhardt and Piwnicki [17] and Lorenci 
and Klippert [18], for example, have discussed the case of electromagnetic propagation in 
moving media. In order to achieve four dimensional geodesic flows, it has been necessary 
to introduce a motion of the medium. 2 There is considerable interest in extending these 
analog models for the kinematical aspects of gravity to include dynamical aspects, i.e. 
considering gravity as an emergent phenomenon [16]. 

The manifestly covariant classical and quantum mechanics introduced by 
Stueckelberg [6] in 1941 has the structure of the Hamilton dynamics with the Euclidean 
three dimensional space replaced by four-dimensional Minkowski space (since all four of the 
components of energy-momentum are kinematically independent, the theory is intrinsically 
"off-shell"). This theory leads to five dimensional wave equations for the associated gauge 
fields (in fact, for wave phenomena in general, such as acoustic or hydrodynamic waves), 
which in the eikonal approximation in the presence of an inhomogeneous medium, provides 
a basis for geodesic motion in four dimensional spacetime. 

This theory has been used to account for the known bound state spectra of the (spin- 
less) two-body problem in potential theories formulated in a manifestly covariant way [21] 
(as well as to study the classical relativistic Kepler orbits [6,22]). As we have remarked 
above, in order for the Stueckelberg-Schrodinger equation of the quantum form of this the- 
ory to be gauge invariant, it is necessary to introduce a fifth gauge field, compensating for 
the derivative of the invariant evolution parameter r [10]. Generalized gauge invariant field 
strengths, f QyP , with q,p = 0, 1, 2, 3, 5 occurring in the Lagrangian to second order generate 
second order field equations analogous to the usual Maxwell equations, with source given 
by the four-vector mattter field current and an additional Lorentz scalar density. Taking 
the Fourier transform of these equations over the invariant parameter, as we demonstrate 
below, one sees that the zero frequency component (zero mode) of the equations coin- 
cides with the standard Maxwell theory. The field equations go over in the same way; the 
equation for the fifth field decouples. The Maxwell theory is therefore properly contained 
in the five-dimensional generalization as the zero frequency component. In the quantum 
case the four-currents are given by bilinears in the wave function containing first deriva- 
tives and the fifth source is the (Lorentz) scalar probability density. The symmetry of the 
homogeneous equations, which can be (9(3,2) or 0(4, 1), depending on the sign chosen 
for raising and lowering the fifth index, is therefore not realized in the inhomogeneous 
equations. Hence, without augmenting the symmetry of the matter fields beyond 0(3, 1), 

2 We remark that Obukhov and Hehl [19] have shown that a conformal class of metrics 
for spacetime can be derived by imposing constrained linear constitutive relations between 
the electromagnetic fields (E,B) and the excitations (D,H), using Urbantke's formulas 
[20], developed to define locally integrable parallel transport orbits in Yang-Mills theories 
(on tangent 2-plane elements on which the Yang-Mills curvature vanishes) . 
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the fifth field, whose source in the Maxwell-like equations is the probability density (or, 
classically,the matter density) in spacetime, can play a special role. There appears to be no 
kinematic basis for choosing one or the other of these signatures; atomic radiative decay, 
for example, contains points in phase space (for radiation of off-shell photons) for either 
type. We note, however, that the homogeneous equations corresponding to the (9(4, 1) 
signature appear, under Fourier transform of the r variable, as Klein-Gordon type wave 
equations with positive mass-squared (physical particles) , but for the 0(3, 2) choice of 
signature, these equations have the "wrong" sign (tachyonic) for interpreting the additive 
term as a mass-squared. We do not find an objective criterion for choosing one or the 
other of these possibilities, and therefore leave open the question of a definitive choice of 
the signature for the five dimensional radiation field at this stage. It is possible that both 
may play a dynamical role. 

The structure of the gauge theory obtained from the non-relativistic Schrodinger equa- 
tion is precisely analogous. The four dimensional gauge invariant field strengths obey 
second order equations, for which the sources are the vector currents and the scalar prob- 
ability density (or, classically, the matter density in three dimensional space). No linear 
transformation can connect the Schrodinger probability density with the vector currents, 
so that the scalar density and the fourth gauge field can play a special role. 

Since the eikonal approximation naturally lowers the dimension of the differential 
equations describing the fields by one, the eikonal approximation to the five-dimensional 
field equations results in four dimensional differential equations. 

In the presence of a non-trivial dielectric structure of the medium, the four dimensional 
field equations resulting from the eikonal approximation can describe geodesic motion in 
four dimensional spacetime without the necessity of adding motion to the medium. We 
emphasize that the underlying manifold, on which the fields are defined, is a flat Cartesian 
space, but that the dynamically induced trajectories move on the geodesies of a pseudo- 
Riemannian manifold. This result forms our basic motivation for studying a generalized 
dynamics of Stueckelberg [6] form in this context. 

In Section 5, we study the eikonal structure of waves in a 5D inhomogeneous medium, 
in which Minkowski spacetime is embedded, both for second order wave equations and a 
wave equation of Schrodinger type, and show that the resulting rays have a precise analog 
with the results of Kline and Kay [13], but in a 4D spacetime manifold with a pseudo- 
Riemannian structure. Kline and Kay [13] show that the rays are geodesic in the metric 
associated with the anisotropic inhomogeneous medium. As for the case of Kline and Kay 
[13], it follows from the existence of a Hamiltonian that the corresponding Lagrangian 
obeys an extremum condition which describes the rays as geodesies. We show that there 
is mass flow along these rays, and that the flow is controlled by generating functions of 
Hamiltonian type, establishing a relation between geodesic flow and a particle mechanics 
of symplectic form. We stress that this construction does not necessary imply that the 
equation is written in a curved spacetime, as for the wave equations of Birrell and Davies 
[23], for example. The quantum theory associated with this equation is defined on a flat 
Minkowski spacetime; the wave functions satisfy Euclidean normalization conditons, and 
four-momentum remains the generator of spacetime translations. It is only in the eikonal 
approximation that the rays emerge as geodesies of a curved space. 



8 



This result is the principal content of this work. 

We close our study here by considering a basis for the structure of the Steuckelberg- 
Schrodinger equation with a second rank Lorentz tensor applied to the kinematic terms. 

The ideas of stochastic processes originated in the second half of the 19th century 
in thermodynamics, through the manifestation of the kinetic theory of gases. In 1905 
A. Einstein [24] in his paper on Brownian motion provided a decisive breakthrough in 
the understanding of the phenomena. Moreover, it was a proof convincing physicists of 
the reality of atoms and molecules, the motivation for Einstein's work. It is interesting 
to note that Einstein predicted the so called Brownian motion of suspended microscopic 
particles not knowing that R. Brown first discovered it in 1827 [24]. The resemblance of 
the Schrodinger equation to the diffusion equation had lead physicists (including Einstein 
and Schrodinger) to attempt to connect quantum mechanics with an underlying stochastic 
process, the Brownian process. 

Nelson [1], in 1966, constructed the Schrodinger equation from an analysis of Brownian 
motion by identifying the forward and backward average velocities of a Brownian particle 
with the real and imaginary parts of a wave function. He pointed out that the basic process 
involved is defined non-relativistically, and can be used if relativistic effects can be "safely" 
neglected. The development of a relativistically covariant formulation of Brownian motion 
could therefore provide some insight into the structure of a relativistic quantum theory. 

Nelson pointed out that the formulation of his stochastic mechanics in the context of 
general relativity is an important open question [1]. The Riemannian metric spaces one 
can achieve, in principle, which arise due to nontrivial correlations between fluctuations 
in spatial directions, could, in the framework of a covariant theory of Brownian motion, 
lead to spacetime pseudo- Riemannian metrics in the structure of diffusion and Schrodinger 
equations. Morato and Viola [25] have recently constructed a relativistic quantum equa- 
tion for a free scalar field. They assumed the existence of a 3D (spatial) diffusion in a 
comoving frame, a non-inertial frame in which the average velocity field of the Brownian 
particle (current velocity) is zero. In this frame the location of the Brownian particle in 
space experiences Brownian fluctuations parametrized by the proper time of the comoving 
observer. They interpreted possible negative 0-component current velocities with what 
they called 'rare events', which are time reversed Brownian processes (a peculiarity arising 
in the relativistic treatment). The equation they achieved this way is approximately the 
Klein-Gordon equation. It is important to note that in the inertial frame they do not 
obtain a normal diffusion. This is due to the fact that their process is stochastic only in 
three degrees of freedom and therefore is not covariant. In this paper we shall study a 
manifestly covariant form of Brownian motion, par In a previous work [27] we introduced 
a new approach to the formulation of relativistic Brownian motion in 1+1 dimensions. The 
process we formulate is a straightforward generalization of the standard one dimensional 
diffusion to 1+1 dimensions (where the actual random process is thought of as a 'diffusion' 
in the time direction as well as in space), in an inertial frame. The equation achieved is 
an exact Klein-Gordon equation. It is a relativistic generalization of Nelson's Brownian 
process, the Newtonian diffusion. In this work we review the relativistic Brownian process 
in 1+1 dimensions [27] where the inclusion of both spacelike and timelike motion for the 
Brownian particle (event) is considered; if the timelike motion is considered as "physical" 
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the "unphysical" spacelike motion is represented (through analytic continuation) by imag- 
inary quantities. We extend the treatment to 3+1 dimensions using appropriate weights 
for the imaginary representations. The complete formalism then can be used to construct 
relativistic general covariant diffusion and Schrodinger equations with pseudo-Riemannian 
metrics which follow from the existence of nontrivial correlations between the coordinate 
random variables. 

Finally, we discuss the possible implications of the process we consider (i.e. a relativis- 
tic stochastic process with Markov property which preserves macroscopic Lorentz covari- 
ance) on the entangled system, where we claim that though fluctuations which exceed the 
velocity of light occur, the macroscopic behavior dictated by the resulting Fokker-Planck 
equation is local. 

II. Equations of Motion and Self-Interaction 

The Stueckelberg-Schrodinger equation which governs the evolution of a quantum 
state over the manifold of spacetime was postulated by Stueckelberg [6] to be, for the free 
particle, 



dr 2M 



■VV (2.1) 



where, on functions of spacetime, p M is represented by —id/dx^ = —id^. 

Taking into account full U(l) gauge invariance, corresponding to the requirement 
that the theory maintain its form under the replacement of ip by e te ° A ip : the Stueckelberg- 
Schrodinger equation (including a compensation field for the r-derivative of A) is [10] 

+ e a 5 (x, T))4> T {x) = * Vv(z), (2-2) 

where the gauge fields satisfy 

a' a = a a + d a A 
under gauge transformations generated by 

for A a differentiable function of {a^, r} = x a ( a = 0, 1, 2, 3, 5), which may depend on r, 
and eo is a coupling constant which we shall see has the dimension l~ x . The corresponding 
classical Hamiltonian then has the form 

K = (P"-eoa^r))^-e a,(x,r)) _ ^ (2g) 

in place of (2.1). Stueckelberg [6] did not take into account this full gauge invariance 
requirement, working in the analog of what is known in the nonrelativistic the 
Hamilton gauge (where the gauge function A is restricted to be independent of time) . The 
equations of motion for the field variables are given (for both the classical and quantum 
theories) by [10] 

Xd a f a (x,r) = e j l3 (x,T), (2.4) 
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where a, (3 = 0, 1, 2, 3, 5, the last corresponding to the r index, j 5 = p is he density of events 
in spacetime, and A, of dimension is a factor on the terms f a ^f a p in the Lagrangian 
associated with (2.2) (with, in addition, degrees of freedom of the fields) required by 
dimensionality. The field strengths are 

fa/3 = ga^ _ ( 2 5 ) 

and the current satisfies the conservation law [10] 

d a j a (x,r) = 0; (2.6) 

integrating over r on (— oo, oo), and assuming that j 5 (x, r) vanishes (pointwise) at |r| — > oo, 
one finds that 

d^(x) = 0, 
where (for some dimensionless rj) [9,11] 

/oo 
drf{x,r). (2.7) 
-oo 

We identify this J^(x) with the Maxwell conserved current. In ref. [28], for example, this 
expression occurs with 

f{x : r) =x»(t)5 4 (x-x(t)), (2.8) 

and r is identified with the proper time of the particle (an identification which can be 
made for the motion of a free particle). The conservation of the integrated current then 
follows from the fact that 

d^f = x»(r)d^\x - x(t)) = ~ x(t)), 

a total derivative; we assume that the world line runs to infinity (at least in the time 
dimension) and therefore the 5-function vanishes at the end points[6,28], in accordance 
with the discussion above. 

As for the Maxwell case, one can write the current formally in five-dimensional form 

j a = x a 5 4 (x(r) -x). (2.9) 

For a = 5, the factor x 5 is unity, and this component therefore represents the event density 
in spacetime. 

Integrating the /(/-components of Eq. (2.4) over r (assuming / m5 (x, r) — > (pointwise) 
for r — > ±oo), we obtain the Maxwell equations with the Maxwell charge e = eo/rj and 
the Maxwell fields given by 



A»(x) = A / a^(x,r)dr. (2.10) 

— oo 
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A Hamiltonian of the form (2.3) without r dependence of the fields, and without the 
terms, as written by Stueckelberg [6], can be recovered in the limit of the zero mode of the 
fields (with a$ = 0) in a physical state for which this limit is a good approximation i.e., 
when the Fourier transform of the fields, defined by 

a M (x,r) = J dsd fl (x 1 s)e- isT 1 (2.11) 

has support only in a small neighborhood As of s = 0. The vector potential then takes 
on the form a^{x,r) ~ Asd M (;r,0) = (As/2n\)A fl (x), and we identify e = (As/27rA)e . 
The zero mode therefore emerges when the inverse correlation length of the field As is 
sufficiently small, and then rj = 2n\/As. We remark that in this limit, the fifth equation 
obtained from (2.4) decouples. The Lorentz force obtained from this Hamiltonian, using 
the Hamilton equations, coincides with the usual Lorentz force, and, as we have seen, the 
generalized Maxwell equation reduce to the usual Maxwell equations. The theory therefore 
contains the usual Maxwell Lorentz theory in the limit of the zero mode; for this reason 
we have called this generalized theory the "pre-Maxwell" theory. 

If such a pre-Maxwell theory really underlies the standard Maxwell theory, then there 
should be some physical mechanism which restricts most observations in the laboratory 
to be close to the zero mode. For example, in a metal there is a frequency, the plasma 
frequency, below which there is no transmission of electromagnetic waves. In this case, if 
the physical universe is imbedded in a medium which does not allow high "frequencies" to 
pass, the pre-Maxwell theory reduces to the Maxwell theory. Some study has been carried 
out, for a quite different purpose (of achieving a form of analog gravity), of the properties 
of the generalized fields in a medium with general dielectric tensor [29]; we discuss this 
study in detail in a later section. Moreover, as we describe below [9] the high level of 
nonlinearity of the theory of the electric charge in interaction with itself may generate an 
approximate effective reduction to Maxwell-Lorentz theory, with high frequency chaotic 
behavior providing the regularization achieved by models of the type discussed by Rohrlich 
[30]. 

We remark that integration over r does not bring the generalized Lorentz force into 
the form of the standard Lorentz force, since it is nonlinear, and a convolution remains. If 
the resulting convolution is trivial, i.e., in the zero mode, the two theories then coincide. 
Hence, we expect to see dynamical effects in the generalized theory which are not present 
in the standard Maxwell-Lorentz theory. In the following, we describe results we have 
obtained in a study of the self-interaction of a classical relativistic charged particle. 

Writing the Hamilton equations 

. fl _dx^_dK . fI _dp» OK 
dr dp^ ' ^ dr dx^ 

for the Hamiltonian (2.3), we find the generalized Lorentz force [9,10] 

Mi^eof^ + fs. (2.13) 

Multiplying this equation by i M , one obtains 

Mi,F = e x M r 5 ; (2.14) 
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this equation therefore does not necessarily lead to a trivial relation between ds and dr. 
The f M 5 term has the effect of moving the particle off-shell (as, in the nonrelativistic case, 
the energy is altered by the electric field). 
Let us now define 

ds 2 

£ = 1+2^ = 1-—, (2.15) 

where ds 2 = dt 2 — dx 2 is the square of the proper time. Since x^ = (p M — eoa M ) /M, if we 
interpret — eoa M )(p M — eoa^) = — m 2 , the gauge invariant particle mass [10], then 

measures the deviation from "mass shell" (on mass shell, ds 2 = dr 2 ). 

We see that the a$ field is strongly associated with the mass distribution; its source is 
the event density (mass density in spacetime), i.e. (in generalized Lorentz gauge d a a a = 0), 

-d a d a a 5 = ef = p. 

We carry out a power series expansion of the Green's function in the neighborhood 
of the parameter r locating the particle on its worldline, and compute the field strengths 
entering into the generalized Lorentz force. This results in the system of equations [9] 
(using a cutoff which we estimate to be of the order of 10~ 23 seconds [31] to avoid 
explicit singularities) 

M(e,e)x» = -- M(£,£ W + F{e)e 2 {x tl + -^—x v x u x»} 

2 1-e 1-e (2 _ lg) 

+ e f ext » v x» + e + 8" „) 

for the orbits in spacetime. We moreover obtain an autonomous equation for the off-shell 
deviation, 

p 

(1 + 2— -) "e - As + Be 2 + Ce - D + Ee 3 + lee = 0. (2.19) 

The coefficients are defined as follows 

2 M 2M(e,e) 4MM)F 5 

F 3 [ 2e 2 2> e 2 F(e) 2e 2 F(e) ' 

2F£, Mv 2F! 2 2 1 ,M . M(e,e) 1 



F 3 2 v 2e 2/ F 3 l-eF 3 y 2e 2 } e 2 F(e)l-e' 

_ 4M(M) 1 , M , „ x 2 2F 1 2 

( - 777777-77(77 + F *) - ^1^3 " (1 _ £)Fs + y/n 



D 



e 2 F{e) 


^3 


4M(e,e) 


Fx 


e 2 F(e) 


F 3 


2F 4 





(2.20) 



E = 4 - 2( ^—^ - —) 
(l-e)F 3 { F 3 2 F 3 > 



F 3 F 3 F 3 2 ' (l-e)F 3 
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We have retained the function M(e, e) in some of the coefficients for convenience, 
although it contains a part linear in e, and would therefore modify the definitions of the 
coefficients somewhat (the form of (4.6) would remain) if we were to redefine them to make 
all the derivative of e terms explicit. 

The basic functions in terms of which the Fi above are defined, extracted directly 
from the Green's function, are 



h(e) 
/s(e) = 



fin 



l+y/e I 
1— v/e I 



+ 



fin 



i+V£ I 

1— v/e I 



£ 2 (l-£) £(l-£) 2 



1 



( £ )§ £ 2 £ 2 (l-£) 



| In 



l + y/i" I 
1— v/e I 



3 
£2 



+ 



e(l-e)" 



(2.21) 



For either sign of e, when £ ~ 0, 



/l(e)-^ + fe+f. 2 + 0(£ 3 ), 
/ 2 (e)--^-^-^ 2 +0(£ 3 ), 
/3(e)~^ + ^+^ 2 + 0(£ 3 ). 



(2.22) 



and therefore the functions defined in (2.21) are smooth near e = 0. We then define the 
auxiliary functions 



1 111 

9i = fi — J2 — 3/3 , <72 = ~ h — 2/3 7 #3 = g/i — 2^ 2 ~~ 2^ 3 

hi = ~~ 2^ 2 ~~ ^ > ^2 = - 2/2 ~ 2^ ' ^3 = (/1 - f'2 - fa) ' 



h 4 = n - n - ri 



(2.23) 



Finally, 



F 2 (e 

Fa(s 
F A (s 
F 5 (e 



gi(e-l) 

g 2 - 2(e - l)h ± 
4^ 

9a + —(e-l)ha 
\h 2 + -(e - l)h 4 
\^-e)ha 



(2.24) 
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and 

F(e) = ^(l-e)+g 3 . (2.25) 

Here, the coefficients of x^ have been grouped formally into a renormalized (off-shell) 
mass term, defined (as done in the standard radiation reaction problem ) as 



where [9] 



M(e, e)=M+^- S) + g 2 ] - e 2 [\f[(l - e) + h 2 ]e, (2.26) 

e 2 = Ze o (2 .27) 
A(2tt) 3 ^ V ; 

can be identified with the Maxwell charge by studying the on-shell limit (compare with 
our discussion above; one concludes that in this case As 2 = \/2ir/j,). If e varies slowly, 
this "renormalized mass" term can have some significance as a measurable mass, but if e 
is rapidly varying, the identification is only formal; its connection with a measurable mass 
is only in the sense of some average over local variations. 

We remark that when one multiplies Eq.(2.18) by x^, it becomes an identity (all of the 

terms except for eof ex t^x l/ may be grouped to be proportional to ^ x **£ + ^ifj)] one must 

use Eq.(2.19) to compute the off-shell mass shift e corresponding to the longitudinal degree 
of freedom in the direction of the four velocity of the particle. Eq.(2.18) determines the 
motion orthogonal to the four velocity. Equations (2.18) and (2.19) are the fundamental 
dynamical equations governing the off-shell orbit. 

It can be shown [9] that Eq.(2.18) reduces to the ordinary (Abraham-Lorentz-Dirac) 
radiation reaction formula for small, slowly changing e and that that no instability, no 
radiation, and no acceleration of the electron occurs when it is precisely on shell. There 
is therefore no "runaway solution" for the exact mass shell limit of this theory. This 
result indicates that in the mass shell limit, the theory is fundamentally different than 
the "standard" theory. The unstable Dirac result is approximate for e close to, but not 
precisely zero. The Dirac instabilities are therefore necessarily associated with at least 
small deviations from mass shell. 

The results of this calculation demonstrate the strong connection between the as field, 
whose source is the mass distribution, and the off-shell mass variations of the particle. In 
the next section we give an explicit geometrical interpretation, constructing a conformal 
metric which replaces the fifth gauge field (as in standard electromagnetism, the scalar 
gauge field cannot be removed in the presence of sources) . 

III. Conformal Equivalence for the Fifth Potential 

In this section, we study the replacement of the fifth potential in a flat space picture 
in terms of a new Hamiltonian which contains only a 4-potential, and takes into account 
the fifth potential with a metric coefficient [11]. The generator of motion is (we assume no 
explicit r dependence in the fields in this section, so that a M (x) is a zero- mode field) 

v nu (Pn - ea^(x))(p v - ea v (x)) 

K r = 9 2M (3 ' 1} 



15 



The conformal structure of g^ v follows from the fact that in (2.3), the contraction of 
Ov ~ ea M (x))(p v -ea v (x)) the Minkowski metric is just K + ea 5 , and hence the conformal 
factor $ is given by he expression preceding (3.5). 

We now verify that this functional of (j> a ,x a ) gives Hamilton equations which are 
equivalent dynamically to those of (2.3). We find 

t =*- = ^ = ^ (3-2) 
where we have defined the auxiliary variable £ by 

Note that the covariant Poisson bracket 

M{£ v ,x x } = {p v ,x x } = 6 v x , 

so that 

{x°,x x } = g° x {i„,x x }=±g° x . 

This type of commutation relation, a generalization [10] of the assumption made originally 
by Feynman, maintains the interpretation of p v as the generator of translations in x x (in 
the corresponding quantum theory, p v — > —id/dx v . We shall return to this point in a more 
general context later. 

This equation gives us the transformation law dx a = g av d£ >v and d£ a = g a ^dx^. The 
second Hamilton equation gives 

dp a _ dK r _ M dg^ ■ ■ ^ da^ (p„ - ea u ) 
dr ~ dx° ~ 2 dx°^ u+g dx° M 1 } 



We now replace x with £ using the transformation law and substitute for p a to obtain 

+ e|S = -^9UA + re&g-t, (3-4) 



The functionals K, K r are different; however, on the physical trajectories they take the 
same numerical value, K. We now show that choosing a conformal metric g^ v = ^{x)rj^ v , 
the Lorentz force derived from K r is the same as the one derived from K. In this case we 
have 

9 "- = ¥R"<"" *W = i + t 1 a .W 

and Eq.(2.8) gives 

Using ^'rf v k, lx ^ tV = ^ we find that Eq. (2.9) is indeed satisfied. This shows that the 
Hamilton equations for the two generators are identical. 
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It is now interesting to examine the geodesic motion of this dynamical system, assum- 
ing the fields are static in r. The Lagrangian in this case is 

M 

L = p^ -H = —g^x»x v + i^a M . (3.6) 
We now make a small variation in x^ 

-> x» + 5x^ 

8S = j dT h^^ x 5x +2g ^ x ^ +ea ^ + d^ x 5x 1 



From the minimal action principal we obtain, by integration by parts of the r derivatives 

jL — ^.m 9 ■ 



(r independence of the field implies 4- = x^-J^) 



multiplying by g Xcr we finally get 

x X = -T^x" + ±i»fl (3.7) 

where /* = g^f afI and = |^ - 

As an example of the application of our result for the conformal metric, let us consider 
the Friedmann-Robertson- Walker universe. 

In the "flat space" Robertson- Walker model (see, e.g., [32]), for the spatial geometry 
characterized by k=0, the metric 

ds 2 = dr 2 - $ 2 (r) (dx 2 + dy 2 + dz 2 ) . (3.8) 
can be brought to the form 

ds 2 = $ 2 (t) (dt 2 - dx 2 - dy 2 - dz 2 ) . (3.9) 
by using the transformation 

r is the time coordinate of a freely-falling object and therefore coincides with our notion of 
universal the r. The function $(r) is often designated by R or a and is the (dimensionless) 
spatial scale of the expanding universe. In the conformal coordinates the time-coordinate 
is therefore related to r, according to the transformation above, through 
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It is interesting to use the Lorentz force in order to achieve the same result. Let us assume 
that a 5 depends on t alone. In this case, the force is 



The relation 



then implies 



i.e., 



* = TT^ (313) 



d<3> ^ ^ 4 e da 5 
dt K dt 



da 5 Kd,h 
"TT 777 v7tT2 / 



dt e v $ 2 

We substitute this in the force equation and multiply by 2t to obtain 



eft 2 K d , 1 x /0 , 

d7 = - 2 M57(^' (314 » 



Finally, putting K = — ^ we arrive at the remarkable result 

dt 1 

which coincides with the transformation (3.11) from the time on the freely falling clock 
r to the redshifted t in the conformal form of the Robertson- Walker metric. We see that 
this t corresponds to the Einstein time satisfying the dynamical Hamilton equations, and 
the conformal factor of the Robertson- Walker metric coincides with the conformal facter 
of the curved space embedding. 

In this construction we have assumed the as field to depend on t alone. The generalized 
Maxwell equations then provide a simple connection between the Robertson- Walker scale 
and the event density. 

The generalized Maxwell equations [10] are 

d a f a = e/, (3.15) 

where j 13 = (j^,p) satisfies dpj@ = d^j^ + d$P = 0, and p is the event density. In the 
generalized Lorentz gauge d a a a =0, we have 

-d a d a a 5 = ef = ep. (3.16) 

Since a 5 depends on t alone (3.16) becomes 

d t V = ep. (3.17) 
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From (3.13), 

a5 = ~^-^ 

so that from (3.17) 

The space-time geometry is related to the density of matter pu through the Einstein 
equations 

= - -g^R = 8ttGT^, (3.19) 

where R^ v is the Ricci tensor, R is the scalar curvature and is the energy-momentum 
tensor. For the perfect fluid model (isotropy implies the T^ v is diagonal) 

= pu^u" +P(g^ + u^u u ). (3.20) 

The (0, 0) component (referring to r) is then 

T TT = 8nGp M . (3.21) 

using the affine connection derived from the metric (3.8) one finds 

G TT = 3^ = 8nGp M (3.22) 

and the (equal) diagonal space-space components are (for example, we write the x,x com- 
ponent 

<?" = "^ 8^ = ^- < 3 - 23 > 

Since T TT is the (0,0) component of a tensor, it follows from (3.11) and (3.21) that 
the matter density in the conformal coordinates is given by 

Pm = ^Pm- (3-24) 

To establish a connection between the density of events in spacetime p and the density of 
matter (particles) p' M , in space at a given time t, we assume that, 



Pm = P At (3-25) 



where At is the time interval (in the conformal coordinates associated with the Stueckelberg 
evolution) in which the events generating the particle world lines are uniformly spread. It 
then follows from (3.25) that 

Pm 
p¥< 



At = ™ (3.26) 
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We now consider two examples. For the static universe, for pm constant, it follows 
from Eq. (3.22) that $ is given by an exponential; it then follows from (3.18) that p is 
constant, so that 

At oc $" 2 . (3.27) 

For the matter dominated universe, where the pressure is negligible [32] , one sees from 
(3.23) that 

2H><I> tt = $ 2 , 

and substituting in (3.18), one finds after changing r derivatives to t derivatives in (3.22) 

that — is constant. It then follows that At oc $ -2 in this well. 

p 

This result implies that, at any given stage of development of the universe, i.e., for 
a given r, the events generating the world lines lie in an interval of the conformal time t 
which becomes smaller as $ becomes large in the order of $~ 2 With the relation (3.11), this 
corresponds, on the other hand, to a narrowing distribution, of order <E> _1 in r, contributing 
to a set of events observed at a given value of the conformal time t. In general, if one 
observes the configuration of a system at a given t, the events detected may have their 
origin at widely different values of the world time r parametrizing the trajectories (world 
lines) of the spacetime events. It would be generally difficult to relate such configurations 
to the configurations in spacetime (at a given r, instead of at a given t) predicted by a 
dynamical theory. However, in this case, we see that the spreading is narrowed for large 
so that the set of events occurring at a given r is essentially the same as the set of 
particles occurring at a given t. The observed configuratoins therefore become very close 
to those predicted by the underlying dynamical model. In the general case, the relation 
between p(r) and pMif) could be very complicated, and it may be difficult to see in the 
observed configurations a simple relation to the dynamical model evolving according to 
the world time. In the static and matter dominated Friedmann-Robertson- Walker model, 
the correspondence between the dynamical theory and observed configurations becomes 
more clear as $ becomes large. 

We have shown that the fifth potential of the generalized Maxwell theory, obtained 
throught the requirement of gauge invariance of the Stueckelberg-Schrodinger equation, 
can be eliminated in the function generating evolution of the classical system by replacing 
the Minkowski metric in the kinetic term by a conformal metric. The Hamilton equations 
resulting from this function coincide with the geodesic associated with this metric, and with 
the Hamilton equations of the original form, i.e., the geodesic equations of the conformal 
metric describe orbits that coincide with solutions of the original Hamilton equations, as 
found in previous work which studied the replacement of an invariant (action-at-a-distance) 
potential by a conformal metric [33]. In this case the geodesic equations are those obtained 
from the conformal geodesic with the addition of a Lorentz force in standard form. 

The Robertson- Walker metric can be put into conformal form. The conformal factor 
of the Robertson- Walker metric can then be put into correspondence with the a 5 field 
of the generalized Maxwell theory and therefore, through its t derivatives (we assume no 
explicit r dependence) with the event density. In both the static and the matter dominated 
models, the set of events generating the world lines of the expanding universe condense 
into progressively thinner slices of the conformal time. 
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We now consider a framework in which one can construct an analog model for gravity 
of much greater generality. 

IV. Eikonal Approximation to Wave Equations 

In this section, we derive the relation between the eikonal approximation to the five 
dimensional generalization of Maxwell theory required by the gauge invariance of Stueck- 
elberg's covariant classical and quantum dynamics to demonstrate, in this approxima- 
tion, the existence of geodesic motion for the flow of mass in a four dimensional pseudo- 
Riemannian manifold. These results provide a foundation for the geometrical optics of the 
five dimensional radiation theory and establish a model in which there is mass flow along 
geodesies. 

We apply this method as well to the interesting case of relativistic quantum theory 
with a general second rank tensor coefficient, as in the conformal case studied above. In 
this case the eikonal approximation to the relativistic quantum mechanical current coin- 
cides with the geodesic flow governed by the pseudo-Riemannian metric obtained from the 
eikonal approximation to solutions of the Stueckelberg-Schrodinger equation. This con- 
struction provides a model in which there is an underlying quantum mechanical structure 
for classical dynamical motion along geodesies on a pseudo-Riemannian manifold. The 
locally symplectic structure which emerges is that of Stueckelberg's covariant mechanics 
on this manifold, and provides the principal result of this work. We apply a technique sim- 
ilar to that used by Kline and Kay [13] in the study of four dimensional (Maxwell) wave 
equations to study the structure of wave equations in five dimensions [34] which follow 
as a consequence of gauge invariance of the covariant classical and quantum mechanics of 
Stueckelberg [6]. Since the eikonal approximation results in the loss of one dimension (a 
high frequency limit), the eikonal method applied to four dimensional wave equations in a 
medium with non-trivial dielectric tensor results in a ray approximation on a manifold in 
three dimensions (Riemannian). In the case of the five dimensional radiation theory, one 
finds a ray approximation corresponding to geodesic flow in a four dimensional pseudo- 
Riemannian manifold. We show that there is a Hamiltonian form for the generation of the 
rays. 

We shall use, in this section, Greek letters for space time indices {ji = 0, 1,2,3) and 
Latin letters to include a fifth index representing the Poincare invariant r parameter in 
addition to the usual 4 spacetime coordinates (e.g., q = 0, 1, 2, 3, 5). The analysis proceeds 
by a generalization of the method [29] of replacing the electric and magnetic vector fields 
by the electromagnetic tensor fields (E 7 B) and the excitation tensor fields (D,H). The 
generalized electromagnetic field tensor is written 

fqi <?2 = dq± a q2 ~ 9 q2 a qi , 

where a q are the so-called pre-Maxwell electromagnetic potentials (as pointed out above, 
the fifth gauge potential as is required for gauge compensation of id$, generating the 
evolution of the Stueckelberg wave function [10]). 
We introduce the dual (third rank) tensor 

h,h hh — c li h h qi <J2 f 

"> — fc J qi <J2 5 
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where e 1 2 3 91 92 is the antisymmetric fifth rank Levi-Civita tensor density. The homoge- 
neous pre-Maxwell equations are then given by 

8 h k hhh = 0, (4.1) 

or, more explicitly = ±d/dr, according to the signature of the r variable, i.e., corre- 
sponding, as we have discussed above, to 0(4, 1) or 0(3, 2) symmetry of the homogeneous 
field equations), 

d 5 e hl25 ^^f qiq2 +d^ l ^^^f qiq2 =0. (4.2) 

We now divide Eq.(4.2) into two cases. In the first, the indices h,fo correspond only to 
space-time indices: 

due"" 5X °fx* + 2cU^ A5 /a5 = - d 5 e^ Xa f Xa + 2d a e^ x f X5 = 0, (4.3) 

where the four dimensional Levi-Civita tensor density This equation, on the 0- 

mode (r independent Fourier components) does not involve any of the usual Maxwell fields 
but only the fifth (Lorentz scalar) electromagnetic field. The second set from Eq.(4.2) 
corresponds to h or/ 2 = 5. It is clear then that all the other 4-remaining indices must be 
space-time indices and we obtain 

e^ Sv d a f Sv = - e^dvU* = 0. (4.4) 

It is this equation that reduces on integration over all r, to the two usual homogeneous 
Maxwell equations. The operation of integrating over all r extracts the zero frequency 
components (in r), which we call the 0-mode. This has the effect of reducing the pre- 
Maxwell system of equations, as we discuss below, to the usual Maxwell equations. With 
appropriate identification of the integrated quantities, the zero mode of the pre-Maxwell 
equations coincides with the Maxwell theory (it is for this reason that the five dimensional 
gauge fields associated with the Stueckelberg theory are called "pre-Maxwell" fields). 

We now turn to the current dependent pre-Maxwell equations. These can be written 

as 

d h n hh =-j l \ (4.5) 

where n^ 2 are the matter induced (excitation) fields (corresponding to H, D in the 4D 
theory). We remark that, restricting our attention to the spacetime components of Eq. 
(4.5), which then reads 

d 5 n» 5 + d„n^ = (4.6) 

we may extract the 0-mode by integrating over all r. Since j k satisfies the five dimensional 
conservation law dkj k = 0, its integral over r (assuming 12 j 5 — > for r — > ±oo) reduces to 
the four dimensional conservation law d^J^ = 0, where J M = j j tJ, (x,r)dT is the 0-mode 
part of j M (this formula for the conserved is given in Jackson [28]). The first term of 
the left side of (4.6) vanishes (assuming that n M5 — > for r — > ±oo), and one obtains the 
form 
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where we may identify the zero mode fields F^, as above, with the Maxwell fields H, D. 
With the zero mode of (4.4), we see that the Maxwell theory is properly contained in the 
five dimensional generalization we are studying here. 

We assume the existence of linear constitutive equations in the dynamical structure 
of the 5D fields in a medium which connects the n tensor-field to the k tensor-fields using 
a fifth rank tensor S which is a generalization of the fourth rank covariant permeabilty- 
dielectric tensor [29] which relates the E, B fields to the excitation fields D, H in the usual 
Maxwell electrodynamics. The constitutive equations have the form 

n hh =8 hi, q , q , qakqiq2q ^ (47) 

antisymmetric in as well as qiq2Q3 (the indices of k have been lowered with the 
Minkowski metric tensor; we shall treat other tensors in the same way in the follow- 
ing). It is useful at this point to distinguish between the space-time elements / Ml/ and the 
elements / M 5, and assume that the tensor introduced in Eq. (4.7) does not mix these fields 
(for n M5 , if £^5 qiq2 q 3 has gig2<?3 = a/?7, then the components of k that enter are of the 
form /Ca/37 = Zafi-ynbf^ 5 only; similarly, for n^, only the components £ M z, a /?5 can occur, 
and k a @ 5 connects only to the components f\ a of the field tensor). As we have pointed 
out above, the vector / m5 is physically distinguished from the antisymmetric tensor f^ u in 
the inhomogensous field equations, since the source terms break the higher symmetry of 
the homogeneous field equations. The assumption that the constitutive equations do not 
couple these components results in a simpler system to analyze, although (as for Hall type 
effects in the non-relativistic theory) it is conceivable that the more general case could 
occur. 

We introduce the new set of fields: 

= \e» vXa fx^ (4.8) 

so that 

n \a = 2E x<ja ^b aP . (4.9) 

On the zero mode, the fields b^ v correspond to the dual Maxwell fields; in this theory they 
play a role analogous to the B fields in the Maxwell theory. In a similar way, the fields 
are analogous to E. The part of the tensor £i 1 1 2 qi q2 qs connecting the ju5 fields is discussed 
below. 

Working with these fields enables us to construct the equations in a form which, as 
we shall show, generalizes the Maxwell theory to a form where the invariant time r plays 
the role of t and spacetime plays the role of space. This analogy helps to interpret the 
physics and it distinguishes between the familiar physical quantites and the new fields 
/ m5 . Substituting these fields in (4.3) and (4.4), we find 

d 5 b^ + d a ev x fx5 = 0, (4.10) 

and 

d a \f a = 0. (4.11) 
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For the spacetime excitation fields we define 



and we get from (4.5), for I2 = 



a 5n M5 _ ^xu d(rhxv = (412) 

where we have used 

e aWS e^ = -2(W£ -6£S%). (4.13) 
To complete the set of equations, we note that for I2 = 5, we get from (4.5) 

dan 5 * = -f. (4.14) 

To obtain a mass-energy conservation law for the fields, we multiply (4.12) by / M 5 and 
(4.10) by h^v, and then combine them, obtaining 

[U5d 5 n^ + h^d T b^} + ^e^dM^hx.) = -ff^- (4-15) 

Assuming the dielectric tensor reduced into the //5 and \iv subspaces is symmetric (the 
relations of n^ 5 to and n^ v to go by the contraction £e; the exclusive property of 
indices of e then imply simple conditions on £ for the symmetry of these forms) we can 
write (4.15) as 

\d 5 [f^ 5 + 1 V^l + \e^d a {f^h Xv ) = -j^fa. (4.16) 
From the Stueckelberg Hamiltonian [10] 

K = 2M ^ ~ eoa ")^ " e ° a ^ ~ e ° a5 ' ( 4 - 17 ) 
it follows from the relativistic Lorentz force [9,10] 

Mx» = e r„x» + e r 5 , (4.18) 

as discussed above, that 

Mx^ = M-^(x^) = x^f^ 5 (4.19) 
ar 

Since, in the Stueckelberg theory, the Hamilton equations imply that 

^ = ^(^-e <), 
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and hence the r derivative in the central equality of (4.19) corresponds to the change in 
the mass-squared (p^ — eoa M )(p^ — e^a^) of the particle. It therefore follows that j^f^s, 
is the rate of mass change of the system. We therefore identify s a = \s ayjXu f^h\ v as the 
analogue of the Maxwell Poynting vector. This Poynting 4-vector is the mass radiation of 
the field. We see, furthermore, that \ [f^n^ 5 + \h^b^ ul \ is the scalar mass density of the 
field (its four integral is the dynamical generator of evolution of the non-interacting field 
[10]). 

We now introduce the eikonal approximation, i.e., set 

fhh( x i r ) = flil 2 ( x ) e *P iK ( T ~ ®( x )) 

for large k. In the absence of sources the 5D-Maxwell equations (4.10), (4.11), (4.12), (4.14) 
take the form (for large k) 

^-£^ ct Wa5 = 0, (4.20) 
Pab^ = 0, (4.21) 

nM 5 + l £ ^p a h Xv = 0, (4.22) 
P«n 5 ° = 0, (4.23) 

where p a = d a ^. 

We now relate the direction of p M to the polarization of the fields. We write the "cross 
product" of n and b (analogous to the cross product of D and B in Maxwell's theory): 

e^xn^b° x = 2n^U Pfll (4.24) 

or 

where we have used (4.20) and (4.23). It is clear that since p^ and the Poynting four- vector 
are cross products of tensors which are not necessarily aligned in the same four-directions, 
they are in general not parallel to each other (in space-time) due to the anisotropy of the 
medium, i.e., the wave normal and radiation flow directions are not, in general, the same. 

The relations (4.20) — (4.23), along with the constitutive relations relating n^ v ,f a \, 
and n^ 5 , f a5: provide relations analogous to those of the standard Maxwell theory char- 
acterizing the possible field strengths of the eikonal approximation in terms of prop- 
erties of the medium. We shall not treat these relations here, but discuss the mass- 
radiation flows, along the rays, on spacetime geodesies in the interesting special case 
where h^ u = b^ u , which is analogous to the case of materials with \i = 1 in Maxwell's 
electromagnetism. This case is interesting since, although the space is empty in the usual 
sense(i.e. E = D , B = H), the dielectric effect involving the / M 5 components can drive the 
radiation on curved trajectories, i.e., the corresponding spacetime can have a non-trivial 
metric structure. 

We multiply (37) by \e al3 ^p p . We then use (39) and (30) to obtain 

nj - PapPfps + Pf3P P f a5 = (4.25) 
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Defining the reduced dielectric tensor Ej 3 as the part of the general dielectric tensor 
which connects only the ah components of the fields, i.e., 



n a 5 



= Sjfps, (4-26) 



the condition (40) then implies that £^fp5 cannot be in the direction of p a (unless it is 
lightlike). We obtain from Eq. (42) 

(S a p - PaP a Sj + Va pP) fa = 0, (4.27) 

where we have chosen the negative sign for the signature of the fifth index, n a ^ = — n Q 5 
[with this choice the flat space limit, for which £ a 13 = S a /3 , Eq. (4.27), with (4.23), 
admits only spacelike p a ; for positive signature of the fifth index, in this limit, p a would 
be timelike]. 

Eq. (4.24) has a solution only if the determinant of the coefficients vanishes (a similar 
calculation in which the field strengths f^ v enter in place of results in the same condition 
on these coefficients, as it must). It is somewhat simpler to work with the eigenvalue 
equation (4.27). Assuming as before that this dielectric tensor is symmetric, we can work 
in a Lorentz frame in which it is diagonal. In this frame we have (for the transformed 
fields) 

Note that in the isotropic case for which all of the £ a are equal, one obtains ppf 135 = 0, 
and the metric becomes conformal, i.e., one obtains the condition 

S- X rf v V^u = -1, 

where r\^ v is the flat space Minkowski metric (—1, 1, 1, 1). 

Multiplying the equation (4.28) on both sides by p a , and summing over a, one obtains 
the condition (p 2 = p M p M ), 

= K = -A- 2 +-A^ + -A^--A- 2 +l. (4.29) 
ti -p z £ 2 - p z £3 ~P £0 - V 

This condition determines, in this case, the Fresnel surface of the wave fronts. 
It then follows that 

dK _ 2p^ +2p » dK (4 30) 

dp^ £> — p 2 dp 2 ' 

Calculating the scalar product of (4.28) and (4.30) one then obtains 
dK 



= -2h r 5 u V {pl)2 ^ —i = o 

KPvJ \J^ 23 (£ l -p 2 ) 2 (£°-p 2 ) 2 dp 2 



(4.31) 
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Multiplying the expression (4.20) for b^ v {h^ v ) by (4.30), the contribution of the second 
term of (4.30) vanishes since the Levi-Civita tensor is antisymmetric; the first term, ac- 
cording to (4.28), is proportional to / a5 , and vanishes for the same reason. It therefore 
follows that 

OK 

t—W = 0. (4.32) 

Since the scalar product of J^r with both and is zero, it is proportional to their 
"cross product" i.e., it is parallel to the Poynting vector. To make the proof explicit, it is 
convenient to define = J^-, Hi = —hoj, Fj = f i5 , and D l = e % ^ k hjk (the space index 
may be raised or lowered without changing sign in our Minkowski metric). In this case, 
the conditions V^h^ v = and V^f^ become 

V°H - V x D = 

-V°f 05 + V • f = (4.33) 
V H = 0, 

where we have used boldface to represent the space components of the vector. In these 
terms, the Poynting vector is given by 



S° = D ■ f 
S = / 05 D + f x H 



05^ . , „ ( 4 - 34 ) 



Taking the cross product of f with the first of (4.33), one obtains 

V°(f x H) = V(f • D) - D(f • V). 

For V° 7^ 0, one may substitute this into the second of (4.34). The / 05 D term, with the 
help of the second of (4.33), cancels, and we are left with 

c0 

S = — V. 

v° 

It then follows that = yrrV^ <, i.e., V^ 1 is proportional to the Poynting vector. For 
the case V° = 0, the second and third of (4.33) imply that 



V f = V H = 0, 



i.e., if V 7^ (the case = is exceptional in the eikonal approximation), it must be 
proportional to f x H. From the first of (4.33), we see that V x D = 0, and if D ^ 0, 
it must be proportional to V. The space part of S^, from the second of (51) is then 
proportional to V. Under these conditions, the time part of vanishes, and therefore we 
again obtain the result that is proportional to S^. If D = 0, then S° = and, since V 
is proportional to f x H, it again follows that is proportional to >S M . 
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From this point, one may follow the same procedure used in the case of Maxwell's elec- 
tromagnetism [13] to obtain the Hamiltonian flow corresponding to the admissible modes. 
The Lagrangian associated with the Hamiltonian (4.29) satisfies a minimal principle, from 
which it follows that the Hamiltonian flow is geodesic on this manifold. Replacing in 
(4.27) by (^" 1 ) /3 7 n 75 , one obtains 

{<T 7 - M^pV'K 5 = 0, (4.35) 

where 

M a ^ = S^^r 7 - ^ v (S~ l r 7 - (4.36) 

The condition (4.23) implies that the solutions can lie only in the hyperplane orthogonal to 
p a . The projection of the matrix M^p lx p v (in the indices a, 7) into the three dimensional 
subspace orthogonal to p a is symmetric, and can therefore be diagonalized by an orthogonal 
(or pseudo-orthogonal) transformation in three dimensions. In fact, the Gauss law and a 
gauge condition restrict the polarization degrees of freedom to three [10] (the eikonal 
approximation is far from the zero mode, which corresponds to the Maxwell limit, for 
which only two polarizations survive) and hence one finds three geodesies. 

For p a timelike, one can choose a Lorentz frame in which the eigenvalue condition has 
the form 

{5% - M- (p' ) 2 K 5 = 0, (4.37) 

and for p a spacelike, 

{8 a 7 - M^G?" 3 ) 2 }™" 75 = 0. (4.38) 

In each of these cases, the matrix can be diagonalized under the little group acting in the 
space orthogonal to p a (leaving it invariant) . For the lightlike case, up to a rotation, p a has 
the form (p, 0, 0, p) . The remaining matrix may then be diagonalized under SO (2) rotations, 
to obtain just two geodesies, corresponding to the polarization states of a massless Maxwell- 
like theory. This special limiting case will be investigated in detail elsewhere. 

With the same procedure as applied to the Maxwell case [13] one finds the symplectic 
structure of the flow of matter in space time. 

It has been shown by Kline and Kay [13], that for the three dimensional Maxwell 
case, the Hamilton equations resulting from the eikonal coincide with the geodesic flow 
generated by the resulting metric (recall that the direction of momentum associated with 
both eigenstates is the same); a similar proof can be applied to the 4D cases we have 
studied here. 

The study of the five dimensional wave equations we have carried out above provide 
an interesting example of the construction of an analog gravity. We have emphasized the 
important special case where the dielectric tensor in the Maxwell part of the constitu- 
tive equations is trivial, and sufficient curvature is introduced through coupling to the 
fifth component alone. We now turn to a study of a quantum equation which results in 
gravitational physics in the ray approximation, providing a quantum theory which may 
underlie the observed classical gravitational fields. In a later section, we give a mecha- 
nism, following Nelsonfl], based on correlations in relativistic Brownian motion, providing 
an interpretation for the structure of this equation. 
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One could think of applying the eikonal method to a Schrodinger equation in a medium 
which is not isotropic, for example, in a crystal with shear forces [35] , with locally varying 
band structure (as in a crystal under nonuniform stress, or near the boundaries or impuri- 
ties); the second order derivatives in the Schrodinger equation then appear multiplied by 
a "mass matrix." 

The rays are directly associated with the (probability) flow of particles. The eikonal 
eigenvalue condition is one dimensional in this case, since the field is scalar. For an 
analog of this structure (corresponding, for example, to a distribution of events periodic in 
both space and time) in four dimensions described by a relativistically covariant equation 
of Stueckelberg-Schrddinger type, the metric one obtains is a spacetime metric, and the 
geodesic flow is that of the quantum probability for the spacetime events (matter) described 
by the Stueckelberg wave function. We propose a study of the equation 



* JVOr) = -^d»g^{x)d^ T {x), (4.39) 



where g^x) is assumed to be symmetric, and is somewhat analogous to a gauge field. We 
shall refer to it as a tensor gauge field. We assume no explicity r dependence in g fIU (x) 
in this work. The Schrodinger current (satisfying the five dimensional conservation law 
(2.6) )is then 

3r{x) v = -^(Vv*<7m^Vv - W<7^VV*). (4.40) 

In the eikonal approximation, for which the frequency associated with r (essentially the 
total mass of the system [21]) is large, we assume a form for the solution 



ip ~ e 



— i(KT— t/kS) 



where S is the eikonal phase. 

One obtains the condition, for k 



^iTPvPv = h (4-41) 

where 

Pfi = d^S 

analogous to the Fresnel surface condition (4.29) for the optical case. 
We define 

2M y 

as the generator of motion in the corresponding classical dynamics. It is clear that dK/dp^ 
is in the direction of j£. This implies that K is the operator of evolution for the dynamical 
flow of particles, corresponding to the rays. It follows from the Hamilton equations that 
the flow is geodesic, where g^ v is the metric for this manifold. 
As in the discussion of the conformal case, we see that 



x^ = 



OK = g^p v 
dp^ M ' 
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so that 



It then follows that 

K = ^Mx x x a g Xa . 

If we assume that in the flat asymptotic limit the particle is on-shell, since K conserved, 
with the value K = — 4^ everywhere, we would obtain 

dr 2 = —g\ a dx x dx a . 

If we identify the world time interval dr with a measure of length (as in free fall), then we 
can understand g^ v as a metric on a psuedo-Riemannian manifold. As for the conformal 
case, however, the function p u induces canonical translations on x v . 

We emphasize that (4.39) is an equation on flat Minkowski space, and the tensor 
gauge function g^ v is given as a second rank tensor under the Lorentz group. It is only 
in the eikonal approximation that one finds geodesic flow governed by g^ v and its affine 
associated connection with which the manifold is naturally lifted to a curved space, and 
the structure may become covariant under local diffeomorphisms (general covariance). The 
equivalence principle then appears as a characterization of the tangent space of this new 
manifold. 

With this equation, one can study a quantum theory which underlies classical gravity 
and, in particular study the quantum behavior in the neighborhood of singularities of the 
metric for which smooth eikonals may not exist. In the case that the eikonal approximation 
is valid, the ray approximation provides the geodesies of the corresponding gravitational 
field, and assures that the probability flow is along the geodesic. 

V. An Interpretation Based on Relativistic Brownian Motion 

In this section, we shall discuss relativistic Brownian motion, and follow Nelson for 
the construction of the Stueckelberg-Schrodinger equation with tensor gauge coupling. In 
this framework, the tensor coupling arises from correlations between spacetime dimensions 
in the underlying Brownian processes. In the following, we pose and solve some of the 
difficulties in achieving a definition of relativistic Brownian motion. 

Brownian motion, thought of as a series of "jumps" of a particle along its path, neces- 
sarily involves an ordered sequence. In the nonrelativistic theory, this ordering is naturally 
provided by the Newtonian time parameter. In a relativistic framework, the Einstein time 
t does not provide a suitable parameter. If we contemplate jumps in spacetime, to ac- 
commodate a covariant formulation, a possible spacelike interval between two jumps may 
appear in two orderings in different Lorentz frames. The introduction of proper time as 
a parameter for the relativistic Brownian process (RBP) is not possible since the second 
order correlations in the simplest case (i.e. for an isotropic homogeneous process with a 
diffusion constant a 2 ) have the form, for each /j,, 

E(Ax^Ax^) = 2a 2 As (5.1) 
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for each /i; however, summing over /i, 



E(Ax ll Ax^) = As 2 oc As, (5.2) 

where the first equality is by the definition of proper time and the second equality is due 
to the Brownian property expressed in Eq.(5.1). There is an obvious contradiction. We 
therefore adopt the invariant parameter r as the dynamical variable for the Brownian 
process, first suggested by Stueckelberg [6]. 

The interpretation of an event going backwards in time as the antiparticle was given 
first by Stueckelberg and was later used by Feynman; it is now an accepted concept. 
In Feynman's perturbative formulation of quantum electrodynamics [6], pair annihilation 
and creation occurs at points which are sharp vertices at the transitions. However, on a 
smooth worldline describing the pair annihilation process, as described by Stueckelberg 
[6], there are segments in which the event goes faster than light speed (either forward 
or backward in t). In the formulation of our relativistic Brownian process, such sectors 
appear to play an important role; the occurrence of such states of motion is dictated by 
the demand of achieving a Lorentz invariant operator (more explicitly, the d'Alembertian) 
in the relativistic diffusion equation. 

A second fundamental difficulty in formulating a covariant theory of Brownian motion 
lies in the form of the correlation function of the random variables of spacetime. The 
correlation function for the usual isotropic non-relativistic Wiener (.e.g. [1]) process is 
given by 

dw l (t)dw J (t') = a 2 5 %J 5(t-t'), 

for i,j = 1,2,3, where dw l corresponds to the Brownian random part of the Langevin 
evolution (f3 l corresponds to a smooth drift) 

dx l = Fdt + dw l (5.3) 

A straightforward covariant generalization to the relativistic case is 

dw^(T)dw v (r') = a 2 r]^ v d(T - r') (5.4) 

where rf v = diag(— 1, 1, 1, 1) is the Minkowski metric. It then follows that dwo(r)dw (T) < 
0, which is impossible. Let us consider, however, a process which is physically restricted 
to only to spacelike or timelike jumps. One may argue that Brownian motion in spacetime 
should be a generalization of the non-relativistic problem, constructed by observing the 
non-relativistic process from a moving frame according to the transformation laws of special 
relativity (one could, alternatively, argue that as an idealization of a collision model, the 
"jumps" should be timelike; we shall consider both possibilities in the following). Hence 
the process taking place in space in the non-relativistic theory would be replaced by a 
spacetime process in which the Brownian jumps are spacelike. The pure time (negative) 
self-correlation does therefore not occur. In order to meet this requirement, we shall use a 
coordinatization in terms of generalized polar coordinates which assures that all jumps are 
spacelike. A corresponding distribution for such a relativistic Brownian probability density 
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could be, for example, of the form e «^ , where \i is the invariant spacelike interval of the 
jump. This is a straightforward generalization of the standard Brownian process in 3D, 

-r 2 

which is generated by a probability density of the form e sa* , where r is the rotation 
invariant (i.e. the vector length) and a is proportional to the diffusion constant. We shall 
refer to this function as the relativistic Gaussian. 

As we shall see, a Brownian motion based on purely spacelike jumps does not, however, 
yield the correct form for an invariant diffusion process. We must therefore consider the 
possibility as well that, in the framework of relativistic dynamics, there are timelike jumps. 

CT 2 

The corresponding distribution would be expected to be of the form e~ ra? , where a is the 
invariant interval for the timelike jumps, and b is some constant. By suitably weighting 
the occurrence of the spacelike process (which we take for the present discussion to be 
"physical", since its nonrelativistic limit coincides with the usual Brownian motion) and 
an analytic continuation of the timelike process, we show that one indeed obtains a Lorentz 
invariant Fokker-Planck equation in which the d'Alembert operator appears in place of the 
Laplace operator of the 3D Fokker-Planck equation 3 . One may, alternatively, consider the 
timelike process as "physical" (as might emerge from a microscopic model with scattering) 
and analytically continue the spacelike ( "unphysical" ) process to achieve a d'Alembert 
operator with opposite sign. 



5.1 Brownian motion in 1+1 dimensions 



We consider a Brownian path in 1 + 1 dimensions generated by a stochastic differential 
(analogue to the Langevin equation and Smoluchowsky process [36]), of the form 



dx»{r) = pi*(x(T))dT + dw»(r), (5.5) 

where dw is a random process which is a relativistic generalization of the Wiener process, 
and /3 M is a smooth deterministic field (the drift). 

We start by considering the second order term in the series expansion of a function of 
position of the particle on the world line, /(x M (r) + Aa; M ), involving the operator 

= Ax^Ax v ^-^-. (5.6) 
dx» dx u 

We have remarked that one of the difficulties in describing Brownian motion in spacetime 
is the possible occurrence of a negative value for the second moment of some component 
of the Lorentz four vector random variable. If the Brownian jump is timelike, or spacelike, 
however, the components of the four vector are not independent, but must satisfy the time- 
like or spacelike constraint. Such constraints can be realized by using parameterizations 
for the jumps in which they are restricted geometrically to be timelike or spacelike. We 



3 The d'Alembert operator alone does not, as an evolution kernel, assure that an initial 
positive density remains positive; additional conditions must be imposed. We thank Phillip 
Pearle for a discussion of this point. 
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now separate the random jumps into space-like jumps and time-like jumps accordingly, 
i.e., for the spacelike jumps, 



Aw 1 = i/ucosha , Aw° = jusinha (5-7) 

and for the timelike jumps, 

Aw 1 = asinha , Aw° = iacosha (5-8) 

Here we assume that the two sectors have the same distribution on the hyperbolic variable. 
We furthermore assume that fi,a are generated by a relativistic Gaussian distribution, 
working in a Lorentz frame where the a distribution is assumed to be independent of /i, a 
and is uniformly distributed on the restricted interval [— L, L] (see discussion below) where 
L is arbitrary large. Therefore, in this frame < Aw^ > is zero (this is true in all frames; 
see discussion in Section 5) and we pick a normalization such that (for any component) 

< Aw n >oc Ar^ so to first order in At the contribution to < O > comes only from 

< Aw^Aw" >. 

For a particle experiencing space-like jumps only, the operator O takes the following 
form: 

q2 q2 q2 

O spacelike = ^ 2 [ cos h 2 a^2 + 2 sinh a cosh a^-^ + sinh 2 a— ] (5.9) 
If the particle undergoes time-like jumps only the operator O takes the form: 

q2 q2 q2 

Otimehke = cr 2 [sinh 2 a-^ + 2 sinh a cosh + cosh 2 a-^ ] (5.10) 



x 



Since fi, a and a are random processes, the average value of the operator O is the 
sum of the two averages of Eq.(5.9) and Eq.(5.10). A difference between these two aver- 
ages, leading to the d'Alembert operator can only be obtained by considering the analytic 
continuation of the timelike process to the spacelike domain, choosing fx 2 = —a 2 . 

This procedure is analogous to the effect, well-known in relativistic quantum scattering 
theory, of a physical process in the crossed (t) channel on the observed process in the direct 
(s) channel. For example, in the LSZ formulation of relativistic scattering in quantum 
field theory (e.g., [37]), a creation operator in the "in" state may be moved to the left in 
the vacuum expectation value expression for the S'-matrix, and an annihilation operator 
for the "out" state may be moved to the right. The resulting amplitude, identical to the 
original one in value, represents a process that is unphysical; its total "energy" (the sum of 
four-momenta squared) now has the wrong sign. Assuming that the S'-matrix is an analytic 
function, one may then analytically continue the energy-momentum variables to obtain the 
correct sign for the physical process in the new channel. Although we are dealing with an 
apparently classical process, as Nelson has shown, the Brownian motion problem gives rise 
to a Schrodinger equation, and therefore contains properties of the differential equations 
of the quantum theory. We thus see the remarkable fact that one must take into account 
the physical effect of the analytic continuation of processes occurring in a non-physical, in 
this case timelike, domain, on the total observed behavior of the system. 
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In the timelike case, the velocity of the particle Aw 1 / Aw° < 1. We shall here use the 
dynamical association of coordinate increments with energy and momentum 

Aw° Aw 1 

E = M— — p = M—- , 5.11 

At At v ' 

so that 

a 2 =(^r) 2 (E 2 -p 2 ), (5.12) 

where M is a parameter of dimension mass associated with the Brownian particle. It then 
follows that E 2 — p 2 = cr 2 > 0. For the spacelike case, where p/E > 1, we may con- 

sider the transformation to an imaginary representation E — > iE' and p — > ip', for E',p' 
real (this transformation is similar to the continuation p — » ip' in nonrelativistic tunnelling, 
for which the analytic continuation appears as an instanton), but E 2 —p 2 ^ p' 2 — E' 2 > 0. 
In this case, we take the analytic continuation such that the magnitude of a 2 remains 
unchanged, but can be called — \i 2 , so that E' 2 —p' 2 = p 2 with \i imaginary. The spacelike 
contributions are therefore obtained in this mapping by E, p — > iE, ip and a — > ifx, assuring 
the formation of the d'Alembert operator when the timelike and spacelike fluctuations are 
added with equal weight (this equality is consistent with the natural assumption, in this 
case, of an equal distribution between spacelike and timelike contributions) . The preserva- 
tion of the magnitude of the interval reflects the conservation of a mass-like property which 
remains, as an intrinsic property of the particle, for both spacelike and timelike jumps. 
As mentioned before, one recalls the role of analytic continuation in quantum field theory; 
for the well known Wick rotation (e.g., [38]), however, in that case, only the 0-component 
is analytically continued and no clear direct physical idea or quantity is associated with 
it. In the RBP the identification of the imaginary 4-momentum is dynamical in origin. It 
is due to the Lorentz structure of spacetime, which distinguishes the transitions > 1 
from those with A^ < \ Though one may object to the association of Ax 11 with a dy- 
namical momentum (since the instantaneous derivative is not defined for a Brownian 
process) the Brownian motion could be understood as an approximation to a microscopic 
process, just as it appears in Einstein's work in 1905 [24], where it is assumed that the 
Brownian motion is produced by collisions. The effective conservation of E 2 — p 2 as a real 
quantity in both timelike and spacelike processes suggests that it is a physical property 
which preserves its meaning in both sectors. 

With these assumptions, the cross-term in hyperbolic functions cancels in the sum, 
which now takes the form 

° = " 2 4-|ri (5.13) 

Taking into account the drift term in Eq.(6.2.1), one then finds the relativistic Fokker- 
Planck equation 

where d/dx^ operates on both /3 M and p. 
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We see that the procedure we have followed, identifying a 2 = —[i 2 and assuming equal 
weight, permits us to construct the Lorentz invariant d'Alembertian operator, as required 
for obtaining a relativistically covariant diffusion equation. 

To see this process in terms of a higher symmetry, let us define the invariant k 2 = 
E 2 — p 2 > for the timelike case; our requirement is then that E 2 — p 2 = —k 2 for the 
spacelike case. In the framework of a larger group that includes k as part of a three vector 
(E, K,p), the relation for the timelike case can be considered in terms of the invariant of 
the subgroup 0(1, 2), i.e., E 2 — k 2 — p 2 . The change in sign for the spacelike case yields 
the invariant E 2 + k 2 — p 2 ; we designate the corresponding symmetry (keeping the order 
of E and p) as 0(2,1). These two groups may be thought of as subgroups of 0(2,2), 
where there exists a transformation which changes the sign of the metric of the subgroups 
holding the quantity n 2 constant. The kinematic constraints we have imposed correspond 
to setting these invariants to zero (the zero interval in the 2 + 1 and 1 + 2 spaces) . 

The constraint we have placed on the relation of the timelike and spacelike invariants 
derives from the properties of the distribution function and the requirement of obtaining 
the d'Alembert operator, i.e, Lorentz covariance of the diffusion equation. It appears that 
in order for the Brownian motion to result in a covariant diffusion equation, the distribution 
function has a higher symmetry reflecting the necessary constraints. The transformations 
E — > iE' and p — > ip' used above would then correspond to analytic continuations from one 
(subgroup) sector to another. We shall see a similar structure in the 3 + 1 case, where the 
groups involved can be identified with the symmetries of the U(l) gauge fields associated 
with the quantum Stueckelberg-Schrodinger equation. 

We now investigate the Lorentz invariance of the process and the correlation func- 
tions. The averaging operations are summations with weights (probability) assigned to 
each quantity in the sum. The sums in the continuum are, of course, expressed by inte- 
grals. If we wish to assign a relativistic Gaussian distribution function then the hyperbolic 
angle integration is infinite unless we introduce a cutoff. The question then arises whether 
our process is invariant or not. 

We will show that we can use an arbitrary non-invariant (scalar) probability distri- 
bution (for example, a cutoff on the hyperbolic angle) and still obtain Lorentz invariant 
averages, using the imaginary representations of the 'unphysical jumps'. For example, 
< Aw 11 > stands for a summation with a scalar weight (given by the density) over all 
the vectors Aw* 1 , in the domain. It is therefore a vector. Moreover under the imaginary 
representation of spacelike increments relative to the timelike ones (here we assume the 
timelike jumps physical), Aw^ is a simple vector function over all spacetime which has the 
following form: 

/± w v _ / Am/ m Aw ,fl timelike ^ 
\ iAw ,fl Aw'^ spacelike ' 

where the Am/ m are real. The quantity < Aw^ > (formally written as a discrete sum) is 
given therefore by 

<Aw>*> = J2 P ( Aw ) Aw '*= P'i^w^Aw'" + i P'i^w')^' 1 " 

timelike spacelike (5.16) 

= < Aw' 11 >timelike +i < Am/'" > S pacelike 
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where P(Aw) (or P'(Aw')) is the probability(weight) of having the vector Aw (or Aw'). 
The two vectors in the last equality in Eq.(5.16) are just normal Lorentz vectors. If we now 
pick a distribution in a given frame for which the average of each of them (independent 
of the other) is zero then < Aw^ >= is true in all frames since the 0- vector is Lorentz 
invariant. 

Constructing the second correlation, with the assumption of no correlation between 
spacelike jumps and timelike jumps, we find 

< Aw^Aw" > = J2 p ( Aw ) Aw " Aw " = Yl P'(^w')Aw'^Aw" y + 

timelike 

+ i 2 P'{&w')Aw'»Aw' v 

spacelike 

= - a^ v 

timelike spacelike 

oc rf v DAr, 

where a^ v is the correlation tensor in each case. From the definition of Au/ M (a four 
vector) it follows that a^ v are real Lorentz tensors. The last equality in Eq.(5.17) is a 
demand that could be achieved for the general 1 + n case, by assuming that in a given 
frame there is an invariant Gaussian distribution where the distribution is uniform in all 
angles and that there is a cutoff in the hyperbolic angle. The sum of the two covariant 
tensors (each a result of summation on different sectors) is a Lorentz invariant tensor. 
The higher correlation functions do not interest us since they are of higher order in p and 
therefore in At and do not contribute to the Fokker-Plank equation. 

The mapping given in Eq.(5.15) leads necessarily to a deviation from the standard 
mathematical formulation of Brownian motion. There the probability that a particle start- 
ing at x at time r ending at x' at time t' is equal to the probability that the particle starts 
at x at time r passing through any possible intermediate point x" at time r" < r' and 
going from there to the point x' at time r' . This property is expressed in the Chapman- 
Kolmagorov equation (e.g. [36]), 

p(x, r, x', At') = [ p(x, r, x", r")p(x" , r", x', r')d 4 x" , (5.18) 
Jn 

In the relativistic formulation the vector Aw' = x — x' could be a timelike vector therefore 
resulting in a real valued vector Aw according to the mapping in Eq.(5.15) However, the 
two intermediate vectors AwJ = x — x" and Aw' 2 = x" — x could be spacelike, and take the 
event out of the real manifold into a complex valued coordinate. In this case the Chapman- 
Kolmagorov equation does not hold, and the event may be found outside of the real 
manifold. In order to build a consistent process one must adopt the concept of 'Brownian 
jumps' which could be a result for example of a process in which the event (similar to 
Einstein's original construction) undergoes collisions and for each collision, or 'jump' the 
mapping in Eq.(5.15) holds. Therefore at each point in the physical manifold the event 
may take any increment spacelike or timelike (with a possible complex valued contribution 
to the averages). However, although the vector leading from the initial point , say O, to 
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(5.17) 



the end point, A, may be spacelike and therefore be represented as an imaginary vector it 
is understood that the event arrives at the real spacetime point A, never physically leaving 
the real spacetime. This structure separates the two manifolds, spacetime which is real 
and represents the physical coordinates of the event and a complex space representing the 
processes the event undergoes (virtually) going from one point to another. This structure 
differs in that sense from the mathematical formulation due to Wiener and others, but 
still it can be shown that the process is invariant on the average under decomposition 
into shorter time subprocesses . In other words, we consider the event starting at some 
arbitrary point, and going for some time At. We next decompose the time interval into 
M intervals, so that: 

M 
i=l 

We consider then the expression appearing in the Fokker-Plank equation 

M M 

< Ax»Ax u >=< Q^AxfiQ^Ax?) > . 

i=i j=i 

Since we assume that any two non-equal time jumps are not correlated, i.e. < AxiAxj = 
> for i 7^ j, which leaves only the equal time averages in the sum, 

M M 

< Ax»Ax v >= < Axf A< >= a V" Ar * = Ar 

i=i i=i 

where a 2 is the diffusion constant and we used Eq.(5.17) going from the second piece in 
the equality to the third. 

The notion of "jumps" suggests the consideration of discrete processes, which can also 
be formulated within the relativistic framework and leads, under certain assumptions, to 
a covariant Fokker-Plank equation. For example let us assume a physical process in which 
the "jumps" occur in a very ordered way every tj seconds with a very small time spread 
(i.e. a very small probability that a collision occurs within a time different significantly 
from tj). Then, averaging the 'jumps' over a period r >> tj leads to 

< Aw" Aw" >= Na 2 r]^TjNTj < r < (N + 1)tj 

This result is due to the fact that under our assumptions during the time r, iV single 
'jumps' within separation of each other of tj occurred. The average in Eq.(5.17)does not 
change when r changes in less then tj; however if tj is small then one can replace Eq.(5.17) 
with 

Therefore we recover the standard result for Brownian motion. However there is one very 
important difference which is the fact that r can be taken to be finitely small where in the 
standard Brownian process r can be actually taken to zero. This implies that higher order 
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derivative terms enter into the resulting 'diffusion' equation. For example for an isotropic 
homogeneous Gaussian distribution there will be additional even order derivative operators 
beyond the second order (d'Alembert) with coefficients a n Tj^2n — 1) where n is the (even) 
order of the differential operator. Since both tj and a 2 are small these operators could 
be neglected in general, though there might be special configurations in which their effect 
may be significant. In the following we assume that the tj are very small compared with 
the macroscopic scale and that the 'jumps' are practically ordered with zero spread, thus 
the approximation in Eq.(5.17) is valid and no higher order terms are considered. 

5.2 Brownian motion in 3 + 1 dimensions 

In the 3 + 1 case, we again separate the jumps into timelike and spacelike types. The 
spacelike jumps may be parameterized, in a given frame, by 

Aw° = /usinha 

Aw 1 = u cosh a cos 6 sin ■§ 

o (5-19) 
Aw = ix cosh a sin <p sin •# 

Aw 3 = /j cosh a cos $ 

We assume the four variables /U, a, (p are independent random variables. In addition 
we demand in this frame that $ and <f> are uniformly distributed in their ranges (0, tt) and 
(0, 27r), respectively. In this case, we may average over the trigonometric angles, i.e., 1} and 
(j> and find that: 



2 

2 2 2 tij 

< Aw 1 >s ,9 =< Aw 2 >s ,9=< Aw 3 >s ,9= ^--cosh 2 a 



< Aw° 2 >4>,s = fj, 2 smh 2 a 



(5.20) 



We may obtain the averages over the trigonometric angles of the timelike jumps by replac- 
ing everywhere in Eq.(5.20) 



to obtain 



cosh 2 a <-> sinh 2 a , li 2 — > a 2 



2 

< Aw 1 >^,^=< Aw 2 >0,tf =< Aw 3 ><{>,$= — sinh 2 a 



3 (5.21) 
< Aw° >0,i?= or 2 cosh 2 a, 

where a is a real random variable, the invariant timelike interval. Assuming, as in the 1 + 1 
case, that the likelihood of the jumps being in either the spacelike or (virtual) timelike 
phases are equal, and making an analytic continuation for which a 2 — > —A 2 , the total 
average of the operator O, including the contributions of the remaining degrees of freedom 
fx, A and a is 

d 2 

< O > = (< u 2 >< sinh 2 a > - < A 2 >< cosh 2 a >) —^r 

! dt (5.22) 

- (< fx 2 >< cosh 2 a > - < A 2 >< sinh 2 a >) A 
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If we now insist that the operator < O > is invariant under Lorentz transformations (i.e. 
the d'Alembertian) we impose the condition 



2 i 2 2 2 

< ji >< sinh a > — < A >< cosh a > = 

h 2 .2 • U 2 x ( 5 - 23 ) 

— -(< n >< cosh a >- < A >< sum a >) 

Using the fact that 

< cosh 2 a > — < sinh 2 a >= 1 
, and defining 7 =< sinh a >, we find that 

< a 2 >= I±±r < ^ (5 24) 

3+47 r v y 

The Fokker-Planck equation then takes on the same form as in the 1 + 1 case, i.e., the form 
Eq.(5.14). We remark that for the 1 + 1 case, one finds in the corresponding expression 
that the 3 in the denominator is replaced by unity, and the coefficients 4 are replaced by 
2; in this case the requirement reduces to < /j 2 >=< A 2 > and there is no 7 dependence. 
We see that in the limit of a uniform distribution in a, for which 7 — > 00, 

< A 2 >^< /j 2 > . 



In this case, the relativistic generalization of nonrelativistic Gaussian distribution of the 

form e~~dt is e - ^, which is Lorentz invariant. 

The limiting case 7^0 corresponds to a stochastic process in which in the spacelike 
case there are no fluctuations in time, i.e., the process is that of a nonrelativistic Brownian 
motion. For the timelike case (recall that we have assumed the same distribution function 
over the hyperbolic variable) this limit implies that the fluctuations are entirely in the time 
direction. The limit 7 — > 00 is Lorentz invariant, but the limit 7^0 can clearly be true 
only in a particular frame. 

5.3 The Markov Relation and the 4D Gaussian Process 

In developing the previous ideas leading to the formulation of a RBP, we assumed 
that the probability distribution is consistent with the Markov property expressed in the 
Chapman-Kolmagorov equation. However, for the relativistic Gaussian it is not clear 
whether Eq.(5.18) holds. Therefore we now consider an alternative process, using the 
ideas developed above, resulting eventually in the Klein-Gordon equation. Let us consider 
a 2D Gaussian process generated by a distribution of the form : 

P{W > dT) = 2^ eXp CAW °2Ddr AWl ) (5 - 25) 

This distribution corresponds to a Markov process, a standard normalized Wiener process, 
where D is the diffusion constant. We now use the coordinate representation given in 
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Eq.(5.7) and Eq.(5.8) for the timelike and spacelike sectors to transform the distribution 
function in Eq.(5.25) in both sectors to (use /j 2 in both cases): 

1 (5-26) 



2nDdr v 2Ddr 

where timelike 'jumps' are physical and the measure for both sectors is (idfida Then, using 
Eq.(5.15), we get for the combination of the timelike and spacelike contributions (with the 
appropriate sign) of the averages, say, Awq 2 , 



oo />oo 



. — fi 2 cosh 2a 1 

ex P{ TTFTi )d[ida = —Ddr 

2Ddr iv 



(5.27) 



where we integrated over fx first, using 



/•oo p/ n+l \ 



and then integrated over a using 

da 



2 

oo cosh 2a 



(5.29) 



In a similar way one finds that (using Eq.(5.15)) leading to the negative sign) 

< A Wl 2 >= --Ddr (5.30) 

7T 

Since the probability distribution Eq.(6.5.2) is symmetric in Awi in each sector 
< AwqAwi >= as well as the first moments. Therefore we get in this particular frame a 
d'Alembertian. However, following the methods used above, we see that it is an invariant 
result in all Lorentz frames (though in other frames the distribution may not appear to be 
Gaussian). 

Next we consider the application of the 4D form of Eq.(5.25) 

1 -Aw 2 - Awi 2 - Aw 2 2 - Aw 3 2 
P{W > dT) = An 2 D 2 (dr) 2 eXPi 2Ddr } (5 ' 31) 

with measure p?dji cosh a sin QdQdadcf) for the spacelike sector and fj, 3 dfi sinh a sin 

ed6dad(M 

for the timelike sector. 

However, now calculating < A^o 2 > for the timelike case, after averaging over the 
spatial angles 9 and <fi we find, using Eq.(5.19), 

r0 ° r °° 5 -y 2 cosh 2a 
nDHdr) 2 J } _J eXV{ ^2DdT 



< Awq 2 >= — / / u 5 exp( — — ) cosh asinh adada (5.32) 

('/")" ./l) J-, or, j_ 
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and for the spacelike case we get the same result since the spacelike parametrization of 
Awq 2 is proportional to sinh 2 a and the spacelike volume element is proportional to cosh 2 a. 
Therefore if we use Eq.(5.15), adding the contribution of the two sectors one obtains a 
complete cancellation to zero. In order to avoid this we extend Eq.(5.15) to the form 

Aw^ = Aw/ M , Aw'' 1 timelike 

5.33 

Aw^ = i\Aw'^ , Au/ M spacelike 

Before completing the calculation, we discuss the inclusion of the factor A, in Eq.(5.33). 
Let us consider a classical (i.e. non-stochastic) event with a given value m 2 = Aw M Aw M , 
moving in a timelike direction. It then changes its state of motion and starts moving in a 
spacelike direction; according to Eq.(5.33) m 2 changes into A 2 m 2 . Moreover, though the 
event may move according to a Gaussian distribution which makes no distinction between 
timelike and spacelike motions, the outcome of this motion as represented by the Aw, in 
Eq.(5.15); Eq.(5.33) does distinguish the two phases of motion. We shall see that a specific 
value of A is required for the realization of the Fokker-Planck equation. 

The w manifold is complex and it is a function of the motion on the real manifold 
w' . Our macroscopic (physical) equations are written on the real plane of the w manifold. 
One can then visualize the flow of an event in spacetime similar to a motion of a particle 
in a cloud chamber. There as the particle moves the gas condenses, therefore the particle 
leaves a track. The track itself is not the particle but a result of the actual motion of 
the particle and its interaction with the gas in the cloud chamber. The track in the cloud 
chamber is analogous to the complex representation we use for the 'jumps'. 

We calculate first the expectation < Awq 2 > which is the total expectation, summed 
over the timelike and spacelike sectors. Averaging over the spherical angles 6>, <p we get 
using Eq.(5.32) and Eq.(5.33), 

a 2 !"A 2 , , 2 Z" 00 Z" 00 5 / -^ 2 cosh2a x l2 . l2 , , 

< Awq > = D (dr) / / n exp( — ) cosh a sinh ad[ida 

71 J o J-oo 2Ddr 

8(1 - A 2 1 7 f°° cosh 2 a sinh 2 a , 
-Dot I 5 da 



7r J _ OQ cosh 2a 

where Eq.(5.28) was used in the \i integration leading to the last equality in Eq.(5.34). 
Using 

1 1 

cosh 2 a = -(cosh 2a + 1) sinh 2 a = -(cosh 2a — 1) (6.5.12) 

in Eq.(5.34) and integrating over a we get 

< Aw 2 >= ^~ X ^ Ddr^- (5.36) 

7T 2 

where we used 

r°° da 7T 
cosh 2a 2 

- da TT (5 - 37) 



h = 



cosh 3 2a 4 
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We now calculate the expectation of < Awii 2 >. Averaging over the spherical angles 
6>, (p we get, using Eq.(5.32) and Eq.(5.33), 



< Awi 2 > = — D 2 (dr) 2 [ [ /i 5 exp( - C ° S ^ 2a ) (sinh 4 a - A 2 cosh 4 a)d/ida = 
3tv J J_ 00 2Ddr 

8 „ , f°° (sinh 4 a - A 2 cosh 4 a) , 

= —Ddr / - 3 -da 

Sir J_ 00 cosh 2a 

(5.38) 

Using Eqs.(5.35), (5.29)and (5.37), We get after integration over a, 
<AV>=~[(1-A 2 )f -4(1 + A 2 )] 

In order to obtain the d'Alembertian we insist that < Awi 2 >= — < Awq 2 >, which leads 
to 

A 2 = (5.39) 
3tt + 4 v ; 

Finally, substituting (for example) Eq.(5.39) in Eq.(5.36) we find that 

< Aw^Aw" >= rf v AD dr = rj^Ddr (5.40) 
' 3tt + 4 ' v ' 

where D, is the actual effective diffusion constant defined by 

AD 

D = (5.41) 

3tt + 4 y ! 



VI. Discussion and Conclusions 

We have discussed in the first section the properties of a classical relativistic charged 
particle in the framework of a consistent relativistically covariant classical dynamics. Such 
a theory, without a constraint relating the particle variables E and p admits variations of 
the particle mass from the "mass-shell" (for which an interval of proper time is equal to the 
corresponding interval of the universal world-time); this mass variation plays an important 
role in governing the evolution of the system, and may in fact, macroscopically stabilize 
the spacetime orbit in the presence of the highly singular self-interaction. The work of 
Gupta and Padmanabhan shows that the radiation reaction terms in the Lorentz force of 
the self-interacting system have a geometrical interpretation. In view of this result, one is 
motivated to study a possible connection between such dynamical systems and gravity. 

In the next section, we showed that the a*, field, primarily reesponsible for driving 
the particle off-shell, can indeed be absorbed in an effective conformal metric structure for 
the motion on a manifold. We showed that one can compare this form to the Robertson- 
Walker- Friedmann model, and, through the field equations, obtain a relation for the matter 
density and understand how, in such a world, the Einstein time (in conformal coordinates) 
and r dependence can be constrained to be similar. 
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Following this idea further, we find that the eikonal approximation to wave equations 
for the (generalized) propagation of electromagnetic waves in a medium with non-trivial 
dielectric tensor can result in an analog gravity. 

A mathematically simpler system with this property is that of a Stueckelberg- 
Schrodinger equation with a "tensor gauge" coupling. This equation permits one to study 
a quantum theory on a flat space, for which the eikonal approximation yields a system of 
rays which correspond to the flow of probability along the geodesies of a manifold with 
metric determined by the tensor gauge coupling. Such a system can be lifted by general 
covariance to give the complete structure of general relativity. One can therefore study, in 
this way, the properties of a well-defined canonical quantum theory for which the eikonal 
approximation is classical gravity. The quantum properties of such a system in the neigh- 
borhood of singularities, such as the black hole horizon, may yield interesting information 
on questions such as Bekenstein's conjecture about the spectrum of black hole radiation. 

To understand the structure of a Stueckelberg-Schrodinger equationm of this form, we 
appealed to Nelson's construction of the Schroding equation from the process of Brownian 
motion. To do this, it was necessary to establish a relativic form of Brownian motion, and 
we described and solved the difficulties in achieving such a result. 

In particular, we constructed a relativistic generalization of Brownian motion, using 
the invariant world-time, r, to order the Brownian fluctuations, and separated consider- 
ation of spacelike and timelike jumps to avoid the problems of negative second moments 
which might otherwise follow from the Minkowski signature. Associating the Brownian 
fluctuations with an underlying dynamical process, one may think of 7 discussed in the 
3 + 1 case as an order parameter, where the distribution function (over a), associated with 
the velocities, is determined by the temperature of the underlying dynamical system (the 
result for the 1 + 1 case is independent of the distribution on the hyperbolic variable). More 
generally it is suggestive to consider the possible thermodynamical effects of the 'medium' 
generating the relativistic Brownian fluctuations, following similar steps taken by Einstein 
[Einstein] in his famous work and verify whether any physical effect can be predicted. 

At equilibrium, where dp/dr = 0, the resulting diffusion equation turns into a classical 
wave equation which, in the absence of a drift term K^, is the wave equation for a massless 
field. An exponentially decreasing distribution in r of the form exp — kt would correspond 
to a Klein-Gordon equation for a particle in a tachyonic state (mass squared — k), for 
physical spacelike motion and for physical timelike motion to a particle with mass squared 

K. 

Choosing a cutoff in the hyperbolic angle, one finds a covariant moment and, therefore, 
covariant differential operators. However the underlying process is not invariant; thus one 
can think of a special frame in which the hyperangular distribution is uniformly distributed 
around 0. Boosting breaks the symmetry of the hyperangular distribution, but since the 
the averages are tensor quantities the invariance properties are conserved, and therefore 
the Fokker-Plank equation (leading to the quantum equation) is invariant. This property 
is also used to construct the 4D Gaussian process. 

It was shown that a (Euclidian) Gaussian process with an appropriate (weighted) 
complex representation for the timelike and spacelike random motions can be used to 
achieve the covariant quantum equation, with the assurance that it is Markovian (it is 
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a relativistic generalization of the Wiener process). This leads to a 'cloud chamber-like' 
picture in which the event as it evolves leaves a track (carries a real or an imaginary 
phase), which is a representation of the actual motion, distinguishing the timelike and 
spacelike motion. 

In the classical Stueckelberg theory the timelike (forward or backward) propagation 
is associated with the standard particle or antiparticle interpretation, where spacelike 
propagation is needed whenever one discusses classical pair creation or annihilation (with 
continuous passage from forward to backward motion in time). This suggests that the 
spacelike process may be associated with the annihilation and creation of pairs. Moreover, 
though the resulting macroscopic equation(i.e. on the level of the Fokker-Planck equation) 
is local and causal in the spacetime variables, the underlying microscopic process(i.e., on 
the level of the Brownian fluctuations) is not. It is however local and causal in r even at 
the microscopic level. 

Nelson has shown that non-relativistic Brownian motion can be associated with a 
Schrodinger equation. Equipped with the procedures we presented here, constructing rel- 
ativistic Brownian motion, Nelson's methods can be generalized. One then can construct 
relativistic equations of Schrodinger (Schrddinger-Stueckelberg) type. The eigenvalue equa- 
tions for these relativistic forms are also Klein-Gordon type equations. Moreover one can 
also generalize the case where the fluctuations are not correlated in different directions 
into the case where correlations exist, as discussed by Nelson for three dimensional Rie- 
mannian spaces. In this case the resulting equation will be a quantum equation with a 
local tensor coupling; as we have pointed out, the eikonal approximation to the solutions of 
such an equation contains the geodesic motion of classical general relativity. The medium 
supporting the Brownian motion may be identified with an 'ether' for which the problem 
of local Lorentz symmetry is solved. This study opens up several tracks of possible re- 
search. Nelson [1], discussing the E.P.R. system confronted the fact that such a system 
may be either described by a non-local Markov process or a local non-Markov process. 
The Markov process is simple to implement but Nelson was disturbed by the introduction 
of non-local interaction. However, the non-Markov process is very difficult to apply. The 
RPB developed here may bridge the two possibilities since an ordered (causal) Markov 
process in r may appear to be a non-Markovian (or possibly non-local and certainly non 
causal) process in t. For example for the Gaussian process, looking for the probability of 
finding the event changing its spatial position Ax after At has passed, one may integrate 
Eq.(5.25) over all r. This results however, in Ax aj [ _ At a which is not integrable and therefore 
can not be normalized. This however is not surprising, since the probability of finding the 
event in Ax after At is not well defined (there may be several values of Ax for a given At). 
For example the number of particles, as in Stueckelberg's original construction [6] depends 
on the trajectory through which the point is reached. Defining an appropriate one particle 
probability resulting from the initial process occurring in r demands a restriction of the 
sample space before integrating over r i.e., using the conditional probability restricted to 
processes for which the event's t coordinate is monotonic in r (no pairs are created). 

Finally we would like to point out that generating a covariant quantum equation 
through an RBP leads to a possible relation between quantum mechanics and gravitation. 
In the context of this work, the metric of gravity can appear as an anisotropy in the 



44 



correlations that lead to quantum equations for which the ray, or eikonal approximation, 
corresponds to the classical geodesic flow of general relativity. It furthermore appears 
interesting to generalize Einstein's famous work on this process introducing thermodynamic 
concepts to the resulting geometrical structure of the theory. 

It emerges that the spacetime metric associated with the eikonal geodesic flow, in such 
a theory, would have its origin in correlations in the underlying Brownian process. 
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